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Abstract We review basic theoretical concepts in particle acceleration, with particular em- 
phasis on processes likely to occur in regions of magnetic reconnection. Several new de- 
velopments are discussed, including detailed studies of reconnection in three-dimensional 
magnetic field configurations (e.g., current sheets, collapsing traps, separatrix regions) and 
stochastic acceleration in a turbulent environment. Fluid, test-particle, and particle-in-cell 
approaches are used and results compared. While these studies show considerable promise 
in accounting for the various observational manifestations of solar flares, they are limited by 
a number of factors, mostly relating to available computational power. Not the least of these 
issues is the need to explicitly incorporate the electrodynamic feedback of the accelerated 
particles themselves on the environment in which they are accelerated. A brief prognosis for 
future advancement is offered. 
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1 Introduction 

It is by now well established that the primary source of the energy released in a solar flare 
is in stressed, current-carrying ("non-potential") magnetic fields. This magnetic energy is 
released through a magnetic reconnection process, in which magnetic field lines of opposite 
polarity change their topological connectivity and the field relaxes to a state of lower en- 
ergy, thereby releasing energy in the form of accelerated particles, plasma heating and mass 
motion. 

Magnetic reconnection can occur in several geometries. For example, closed magnetic 
loops embedded into the solar photosphere can collide in the upper solar atmosphere, result- 
ing in a localized region of re latively high magnetic shear, as in the classical Sweet-Parker 
model (see, e.g., ISweetll 19691) . Alternatively, open magnetic field lines can collapse onto 
a closed helmet-like magnetic configuration. This leads to a region of high magnetic shear 
near the top of the closed structure - as evidenced by looptop sources of hard X-rays first ob- 
served by Yohkoh jMasuda etaDll994h - and in the energization of successiv e field lines at 
greater and greater heights, which produces both coronal hard X-ray sources jKmcker et al.l 
|2008a) and a concomitant increase in the separation of the associated hard X-ray and Ha 
footpoints ( iKopp & Pneumanlll976h . Magnetic reconnection can also be accompanied by 
dramat ic large-scale changes in the magnetic field configuration, such as in the "breakout" 
model jAntiochos et alj[l999h . in which a filament, previously restrained by the magnetic 
tension of a set of magnetic field lines exterior to it, is released into interplanetary space as a 
result of reconnection of the constraining field lines. Reconnection can also result in the for- 
mation of magnetohydrodynamic (MHD) shocks, which may be responsible for significant 
particle acceleration exterior to the region of magnetic reconnection itself. 

Whatever the geometry, the fun damental process of magnetic reconnection somehow 
transfers a significant fraction (e.g., lMilleretalJll"99^ of the stored magnetic energy into 
accelerated charged particles, both electrons and ions. Some of the high-energy particles 
may be trapped in the cor ona by plasma turbulen ce (e.g., Petrosian & Liu 2004) and/or by 
magnetic mirroring (e.g., Melrose &Whit3 ll979: Leach 1984; Somov & Kosusi 1997), to 
produce hard X-ray (HXR) looptop sources (e.g..iMasuda et al. , .1994.:,Pe j rosian et al...2C)03 : 
ISui & Holma"i]|2003l : Ijiang et al.ll200^ : lBattaglia & Benzll2006l : lLiull2006l) . Some particles 



may escape the acceleration region into interplaneta ry space and can co ntribute to solar 
energetic particle (SEP) events detected at 1 AU (e.g.. lKrucker et alll2007ah . 

Understanding the conversion of magnetic energy into accelerated particles is one of the 
main goals of solar flare research. It also has attendant implications for our understanding of 
particle acceleration in other astrophysical sites and was indeed a key driver for the design 
of Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI), which measures the 
diagnostic high-energy hard X-ray and gamma-ray radiation produced by the accelerated 
particles with unprecedented temporal, spectral and spatial resolution. 

At a fundamental level, of course, charged particles are accelerated by an electric field. 
Differences in the large-scale structure and temporal evolution of the accelerating electric 
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field allows us to discuss particle acceleration in terms of several basic sc enarios (see Sec- 
tion[3ll. These include large-scale sub-Dreicer acceleration (see, for examplelBenka & Holmanj_ 



1994); super-Dreicer acceleration in a relatively localized 2-D current sheet (e.g., Litvinenko & Somovl 
1993aj); and stocha stic acceleration by MHD turbulence generate d in fast plasma reconnec- 



tion outflows fe.g.. lMiller et al 1 19961 : IPrvadko & Petrosianlll997h 

In Section |4l we report on recent work aimed at extending these basic mechanisms to 
more complex models. These include (1) first-order Fermi acceleration, e ither in a collaps- 
ing rnagnetic trap (|^omov & Ko sugi 1997) or in a fast termination shock ( iTsuneta & Naitd 
ll998l : lAurassetalJl2007h. or bv a super-Dreicer electric field in a sim ple 3-D magnetic 
topology tehark ova & Gor dovskvv '' 2005bl : IZharkova & AgapitovllIoO^ : and (2) accelera- 
tion in more el aborate th ree-dimensional magnetic topologies derived from MHD simula- 
tions (Dalla & Brownindl200a i2008). We also discuss som e recent suggestions of accelera- 
tion by large-scale Alfven waves iFletcher & Hudsonl20p^') and acc eleration invol ving mul- 
tiple reconnection regions or in a turbulent environment dVlahos et a l. 2004; Turk mani et al.l 
l2005l ; lDauphin & Vilmejl2007h . 

In the latter portion of Section |4] we also consider both "local" and "global" aspects of 
the acceleration process. At the local level, accurate simulations of particle trajectories are 
needed to correctly infer the general properties of the acceleration. At the "global" level, we 
recognize that, in view of the large number of accelerated particles involved, the traditional 
"test-particle" approach, in which the electromagnetic environment in the acceleration re- 
gion is prescribed, is not applicable (Section [4.71 - a self-consistent approach requires that 
we incorporate the electric and magnetic fields associated with the accelerated particles 
themselves in current sheets with and without tearing instabilities (Sections 14. 8l and l4~9] l. 

In Sections |5] and [6] we relate the results of the recent modeling efforts to flare obser- 
vations, particularly from RHESSI, and we offer a brief prognosis for future models of the 
particle acceleration process. 



2 Review of Pertinent Observations 

2.1 Hard X-ray lightcurves and spectra 
2.1.1 Lightcurves 

Spatially-integrated hard X-ray light curves obtained by RHESSI confirm many of the tem- 
poral features established by previous observations: sharp increases (bursts) of hard X-ray 
intensity over a relatively short (~0.5-5 s) timescale accompa nied by a more slow ly- varying 
component with a timescale up to a few tens of minutes ( Holman et al.ll20l"l[) . This sug- 
gests that electron acceleration is characterized by two fundamentally distinct timescales: 
a very short timescale associated with the actual acceleration of ambient electrons to high, 
bremsstrahlung-emitting energies and a much longer timescale associated with the duration 
of the acceleration process, a process that, as evidenced by steady hard X-ray emission for a 
substantial fraction of an hour or even longer, can continue to accelerate particles over such 
a sustained period. 

At lower photon energies e, where thermal bremsstrahlung dominates the total emis- 
sion, the emission at higher photon energies is weighted more heavily by plasma at high 
temperatures T. 

As a result, the decrease in conductive cooling time with tem perature - t w T^^l^ for 
coUisionally-dominated conductive processes (see ISpitzedl 19621) - results in the emission 
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at higher energies peaking sooner (lAschwanderj|2007h . and hence the light curve peaking 
progressively earlier with increase in photon energy. At higher energies, where non-thermal 
bremsstrahlung dominates, the relative timing of the emission at different energies depends 
on a balance between the reduced "time-of-flight" for higher energy electrons before they 
impact upon the thick target of the low er atmosphere (which tends to advance high energy 
emissions relative to low energy ones - lAschwanden & Schwartzlll996l : lBrown et allll998h 
and on the decrease in collision freque ncy with ener gy (which tends to delay high energy 
emissions relative to low energy ones - lAschwtmden et al..il997.) . 

2.1.2 Photon and electron energy spectra 

Hard X-ray emission during flares typically exhibits a very steep spectrum at lower energies 
e w 10 keV, in dicative of a thermal process. It must be noted, however, that the frequent as- 
sumption (e.g.. lHolman et al. 2003) of an wothermal pl asma is neithe r expected on the basis 
of sim ple physics nor required by the observations (e.g.. lBrowiJl974h . Indeed. lAschwanderj 
ll2007h has shown that the assumption of an isothermal source is inconsistent with the ob- 
servations of temporal variations of the hard X-ray spectrum in certain flares observed by 
RHESSI. Indeed, in most cases analysis of the variation of residuals (observed hard X-ray 
flux vs. flux calculated using an isothermal spectrum) with photon energy shows that that 
an isothermal model does provide an excellent fit to the observed spectrum; generally a 
broad-band hard X-ray spectrometer shows the (unique) hottest component of a multither- 
mal source. Nevertheless RHESSI's superior resolution and good photometric calibration 
have made it possible to deri ve a two-temperatu re source characterization from a single 
bremsstrahlung spectrum (e.g., ICaspi & LinlboiOh m some cases. 

In most powerful flares the photon spectra reveal double power laws below 300 ke V, with 
the spectr al index below ~30-6 keV being smaller (fl atter) than that above ~ 30-60 keV by 
2-4 units jGrigis & BenzlbOOd : iKrucker et al]l2008ah . For the most energ etic flares, there 
can be discernible emiss ion up to a few hundred MeV (for reviews, see IChupp & RvanI 
l2009l ; IVilmeretal.ll201lb . At photon energies ~500-1000 keV, the spectrum tends to flat- 
ten again, a flattening that has been, at le ast partially, attribu ted to the increasing contribu- 
tion of electron-electron bremsstrahlung dKontar et al.ll2007l) . However, the contribution of 
electron-electron bremsstrahlung to this flattening does not exceed ~0.4, while the observed 
spectral hardening, as observed in a few events with Gj?AA^Ar/PHEBU S, can be as high as 
1.5-2 units at energ ies from 600 keV up to several tens of MeV (see, e.g.. lTrottet et al.l 19981 : 
IVilmeretal.ll 1999 ). 

RHESSI observations reveal a relationship between spectral hardness and time. Values 
of the "local" spectral index y(e) = — J log /(e) /d log e are generally low est around the peak 
of the event - a so-called "soft-hard-soft" evolution originally noted by IParks & Wincklerl 
§969). RHESSFs superior spectral resolution shows that nearly every discernible inten- 
sity peak (which can number u p to a dozen or more in a single flare) has a soft-hard-soft 
structure ( iGrigis &Be"n3l2004h . Also the decrease in spectral index before, and increase 
in spectral index after, this ti me does not occur at the same rate for all photon energies e 
llKontar & MacKinno n' 2005); see Figure 12.11 Rather, a flare typically starts with a spec- 
tral index y that is strongly dependent on energy, with the minimum value of y showing a 
tendency to grow with time. 

Non-thermal emission in the corona is identified in the impulsive phase by its softer 
spectrum (Mariska & McTiernan 1999; Petrosian et al. 2002), consistent with the small col- 
umn depth of the coronal part of the source. The absence of a significant amount of (energy- 
dependent) coUisional losses in this relatively thin target should result in a spectrum two 
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Fig. 2.1 Upper Panel: Temporal variation (4-second cadence) of the count rates in seven front RHESSI seg- 
ments for SOL2002-02-26T10:27 (C9.6) (10:27 UT). The vertical lines show five 20-second accumulation 
intervals for spectral analysis. Lower Panel: Temporal variation of the energy-dependent photon spectral in- 
dex y(e) = —d\ogI(e)ld\os ,e. Each line corresponds to one of the time intervals in the upper plot. From 
iKontar & MacKinnoij i200l . 



powe rs steeper than the target-averaged spectrum ( lBrownll97lllHudsorJl972l : lDatlowe & Lid 
Il973h . 

It should be noted that even if the accelerated electrons have a pure power-law energy 
spectrum Fo{Eq) w ^q"^, for which a characteristic energy does not exist, characteristic en- 
ergies associated with either electron transport and/or radiation physics may produce devi- 
ations from the power-law beh avior in the observed spatially-integrated photon spectrum. 
Further, as elaborated upon bv lKontar et"al] ( l201lh . anisotropy in the mean electron dis- 
tribution, combined with the intrinsic directivity of the bremsstrahlung emission process, 
produces an anisotropic distribution of primary photons. Compton backscattering of pho- 
tons from the photosphere ("X-ray albedo;" Langer & Petrosian 1977; Bai & Ramaty 19^7^ 
not only influences the observed photon spe ctrum, but also has a d iagnostic potential for 
determining the electron angular distribution ( iKontar & Brown|[2006l) . 
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Fig. 2.2 RHESSIhard X-ray photon spectrum for SOL2003-10-28T11:10 (X17.2) (11:06:20-11:06:40 UT). 
The photon flux (plus signs) is fit with a thermal bremsstrahlung contribution from an isothermal plasma 
(dashed curve) plus a double power-law non-thermal electron flux component with a low-energy cutoff (dash- 
dotted curve). The best-fit parameters are shown in the plot. 



2.1.3 Electron numbers 

Analysis of RHESSI data ( iBrown et al.ll2007l : IPennis et al.ll2007h has confirmed the impor- 
tant point that a significant fraction of the energy released in flares is deposited into acceler- 
ated particles. In addition, recent calonme trie analysis of the evolution of a radio spectrum 
from a dense flare (Bastian et al. 1|2003) reveals that a substantial fraction of the energy in 
the energetic electrons can be deposited into the coronal loop, rather than into the chromo- 
sphere; it is estimated that the energy deposited in the corona can approach about 30% of 
the estimated magnetic energy released in the whole flare. 

A s has been pointed out by several authors (e.g.. lBrownll97ll : lTandberg-Hanssen & Emsli3 
Il988h . the total injected electron flux depends critically on the value of the low-energy cut- 
off Ec. A brief discussion of this issue is in order. Historically, a low-energy cutoff Ec was 
assumed simply in order to keep the injected power ^'^^EqFq[Eq) cIEq finite (e.g., iBrownl 
ll97lLlHolman et al.ll2003l) . To determine whether or not such a cut-off is actually required 
by observations, it is best to adopt a non-parametric approach to interpreting the photon 
spectrum 7(e) - i.e ., to infer from /(e) what range of mean ele ctron source funct i ons F (E) 
terown et al.ll2003h allow a statistically acceptable fit to /(e). iKasparova et alJ (|2005*), in 
their analysis of SOL2002-08-20T08:25 (M3.4), have shown that a low-energy cutoff, or 
even a gap, in the mean source electron spectrum exists if the observed spectrum is con- 
sidered as primary bremsstrahlung only; however, such cutoffs and/or gaps disappear if an 
albedo correcti on (for an isotropic primary source) is app lied to the observed photon spec- 
trum (see, e.g.. lKontar et al.ir201 ll) . Further. lEmsliel ( |2003[) has shown that, when allowance 
is made for warm-target effects in the electron energy-loss rate, the injected electron energy 
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Jq E()Fq{Eo) dEo corresponding to a pure power-law photon spectrum can be finite even if 
no low-energy cutoff exists. 

Figure |Z21 shows the photon spectrum for the time interval indicated during SOL2003- 
10-28T11:10 (X 17.2). Assum ing that the non-thermal hard X-ray emission is thick-target 
bremsstrahlung (lBrownlll97lh . one can us e a forward-fitting m ethod, including both ther- 
mal and nonthermal fit components, (e.g., iHolman et ai]|2003h to derive the injected non- 
thermal ele ctron acceler ation rate, differential in energy (electrons s^' keV^'). This method 
(see Kontar et al.ll20I ih also determines an observational upper limit to E^ (as the inter- 
section of the thermal and nonthermal fit components), and hence a lower limit on the 
associated electron number and energy fluxes. Integration over energy provides the total 
electron acceleration rate. For this event, one obtains a total non-thermal electron pro- 
duction rate Fg = 1.7 x 10"^^ s^' above E^ = 28 keV, corresponding to an injected power 
Pe = 1.5 X 10^^ erg s^' above the same energy. For a nominal footpoint area S< 10"* cm^, 
such a particle rate corresponds to an electron flux Eg > 10'^ cm^^ s^' and hence, since the 
average electron velocity v w (c/3) w lO'" cm s^', to a beam density «/, > 10^ cm^-'. This 
is a substantial fraction of the ambient density n w lO'" cm^^^. 

The total electric current corresponding to such a particle acceleration rate, if the ac- 
celeration is unidirectional, is / w 3 x lO'^ A, and the associated current density j = I/S w 
0.3 A cm^^. Such high current densities, especially since they are introduced over a timescale 
of seconds, produce unacceptably large inductive electric and magnetic fields, unless (1) 
the acceleration is near-isotropic, as in stochastic acceleration models (Section [3]4]l, (2) the 
source has very fine structure (Holman 1985), or (3) the beam current is effectively neutral- 
ized b^_a_C0;S£atmlretumcu 

(seelKnight & Sturrock'l977':'Emslie'l980':'Brown & BinghanJl984lLarosa & Emsliell989l : 
lyan den Oord 1990; Zharkova & Gordovskvv 2006). 

Finally, the total number of electrons in the acceleration region = nV w lO-'^, a 
number which corresponds to only a few seconds of accele ration at the rates quoted above. 
Hence the issue of replenishment of the acceleration region ( lEmslie & Henou i 19951) has to 
be accounted for by any acceleration and/or transport mechanism. 

2.2 Lightcurves and energy spectra of y-rays 
2.2.7 Gamma-ray lightcurves 

Due to the limited sensitivity of current gamma-ray detectors, information regarding the time 
scales for ion acceleration is less stringent than that for energetic electrons (Section [2Tb . 

The prompt gamma-ray line emission is known to peak roughly sim ultaneously, withi n 
±1 s, of the hard X-ray continuum emission, based on a few events (e.g.. lKane et al.|[T986l) . 
This near-simultaneity of hard X-ray and gamma-ray emissions indicates that both ion ac- 
celeration to a few MeV, and interaction of th e accelerated ions with the ambient atmosphere 
to produce nuclear lines jVilmer et al]|201lh . can occur within a timescale ~1 s. However, 
in a few flares there are time delays between the appearance of hard X-rays and gamma-ray 
lines (see section 2.2.6 in Vilmer et aLlQlX); this is often interpreted in terms of differential 
transport of electrons and ions between the acce leration and emissi on sites. 

RHESSI obseryations (see, e.g., Figure 3.1 in IVilmer et alJ2oTlh have provided unprece- 
dented information on the temporal evolution of "pure" gamma-ray lines (i.e., lines that oc- 
cur in spectral regions that are relatively uncontaminated by electron bremsstrahlung). Such 
observations indicate the continuous production of very energetic ions (> a few hundreds 
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of MeV) for hours after the impulsiv e phase of the flare , and so pose a formidable constraint 
on the ion acceleration mechanism dVilmer et al.ll201 ll) . 

RHESSI observations have also provided im portant information on the temporal evolu- 
tion of prompt gamma-ray lines dLin et alj2003h . Comparison with the temporal evolution of 
the X-ray flux at 150 keV shows that the temporal evolution of hard X-rays and gamma-rays 
are roughly similar, indicating a common origin of accelerated electrons and ions. However, 
for SOL2002-07-23T00:35 (X4.8), there is a small delay of around 10 s between the max- 
ima of the X-ray and gamma-r ay fluxes. This has been in terpreted in terms o f transport of 
energe tic electrons and ions bv lDauphin & Vilmej ( l2007h . and is discussed in lvilmer et al.l 

ton\) . 

2.2.2 Energy spectra and abundances of ions in flares 

As reviewed bv lVilmeretaTltolll) . interactions of energetic ions in the energy range from 
~1 Mev nucleon^' to 100 Mev nucleon^' produce an extensive spectrum of gamma-ray 
line emission. Narrow gamma-ray lines result from the bombardment of the ambient nu- 
clei by accelerated protons and a-particles, while broad lines occur from the inverse reac- 
tions, in which accelerated heavy nuclei collide with ambient hydrogen or helium nuclei. 
De-excitation gamma-ray line spectra provide information on flare ions of energies above 
~2 Mev nucleon^'. The shape of the ion distribution below this energy is unknown, but is 
comparably important for the evaluation of the total ion energy as is the low-energy cutoff 
to the evaluation of the total electron energy content (Section lTTb (see Ernslie et al. 1997|). 

Quantitative gamma-ray line observations have now been obtained for more than twenty 
events. The accelerated ion spectra are found to extend as unbroken power laws down to at 
least 2 Mev nucleo n~ ' , if a reasonable ambient Ne/0 abundance ratio is used (for details, see 
IVilmeretaUbouh . For the 19 flares observed by the Solar Maximum Mission (SMM) prior 
to the launch of RHESSI, the average spectral index was around 4.3; however, the events 
obser ved with RHESSI tend to have much harder slopes dLin et al.ll2003l : IShare & Murphvl 
l2006h . 

Even though most of the energy contained in ions resides in protons and a-particles, 
the abundances of heavier accelerated ions provide crucial constraints for acceleration pro- 
cesses. Information on the abundances of such flare-accelerated ions was deduced for the 
SMM/GRS (Garrmia-Ray Spectrometer) flares as well as for a few flares observed by RHESSI. 
A noticeable enhancement has been observed in the numbers of a-particles, as well as of 
accelerated ^He isotopes, compared to the standard coronal abundances. Furthermore, ac- 
celerated heavy ions, such as Ne, Mg and Fe, are also generally found to be over-abu ndant 
with respect to normal coronal compositions, as pointed out by Ivilmer et alj ( l201lh . This 
poses another constraint on models for ion acceleration. 

2.2.3 Ion numbers 

Since the first detection of solar flare gamma-ray lines in 1972 ( IChuppetalJl973h . informa- 
tion on the energy content in ions has been obtained for more than 20 events. The analysis 
carried out using 19 events observed by SMM/GRS and a few events observed with RHESSI 
shows that the energy contained in >1 MeV ions ranges from 10^^ to 10^^^ erg, comparable 
to the energy contained in sub-relativistic electrons. However, there are still large uncer- 
tainties in the determination of these quantities and a large dispersion of the ratio of electron 
(^■^ 20 keV) to ion (~ 1 Mev nucleon^' ) energy contents from one flare to another (for details, 
see IVilmeret alj|201ll) . 
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The issue of relative importance of electron and ion acceleration can also be addressed 
by comparing the production ratios of ions to quasi-relativis tic electrons (f^SOO keV) . 
Such an analys is has been perfor med using both S MM/GRS ( lForrestlll98a IChupplll984 
IChupp&BenjI 19941) and RHESSI jShih et alj|2 009b') observations. Good correlations have 
been found between the total energy content in protons above 30 MeV and the total en- 
ergy content in electrons above 300 keV, suggesting that high-energy electrons and ions are 
linked through a c ommon acceleration processes. It should be noted, however, that in the 
statistical study bv lShih et a l. (2009a), a correlation between the number of protons with en- 
ergies above 30 MeV, the electron number above 50 keV, and the soft X-ray Geostationary 
Operational Environmental Satellite ( GOES) class exists, but only above a certain threshold 
in the level of ion production. 



2.3 Geometry of hard X-ray and gamma-ray sources 

The imaging capabilities of RHESSI, combined with its high spectral resolution, have al- 
lowed us to resolve in detail coronal sources (e.g.. lKrucker et al. 2008a) and footpoints oc- 
curring in the same flare (see Figure l23l > and hence to study the acceleration processes that 
give rise to such sources. 

The coronal source often appears before the main flare hard X-ray increase and the 
appearance of footpoints. In the impulsive phase, the coronal hard X-r ay emission is, gen- 
erally, well correlated i n both time and spectrum with the footpoints jEmslie et al.ll2003l : 
lBattaglia&Benzll2006l) . These observations suggest strong coupling between the corona 
and chromosphere during flares, a coupling that is presumably related to transport of accel- 
erated particles from one region to the other. 



2.3.1 Differences infootpoint spectral indices 

An interesting RHESSI observation i s the approximate e quality of the spectral indices in 
different footpoints of the same loop. lEmslie et alJ ( l2003l) reported differences ziy w 0.3 — 
0.4 between the spectral indices of the two dominant footpoints in SOL2002-07-23T00:35 
(X4.8). In a few smaller events analyzed by Battaglia & Benz ( 2006), Zi y is even smaller and 
indeed is significant only in one out of five cases. However, other observations of X-class 
flares reveal a much stronger difference (up to 5) between the spectral indices of footpoints 
of the same loop (Takakura et al. 1995; Battaglia & Benz 2006; Krucker et al. 2007b). 

More surprising than the occasional difference between the spectral indices of foot- 
points are the differences (or lack thereof) between the spectra of the coronal source and the 
footpoints. The spectral index difference Ay = 2.0 predicted llBrownlll97lb between thin- 
and ( collision-dominated) thick-target sources is not observed in most cases ( iHolman et al.l 
I20T 1[). Differences 4 y smaller than 2 could in principle be interpreted by invoking a coronal 



target that is intermediate between thin- and thick-target conditions; such spectral differ- 
ences could also be partially accounted for by the transport effects of precipitating elec- 
trons, e.g., by the effect of a self-indu ced electric field dZharkova & Gordovskyy 200^. 

which 



trons, e.g., by tne eiiect or a seli-indu ced electric tield (iZ,narKova & uordovskyy 
However, 4y ranges from 0.59 to 3.68 jGrigis & Benzll2006l ; iKrucker et alj|2007bh . 



includes values beyond the capability of simple models to explain. For example. Figure |Z4| 
(after Battaglia & Benz 2006) displays the photon energy spectra of two footpoints and a 
coronal source, obtained for the same time interval during the flare SOL2005-07-13T14:49 
(M5.0). The footpoint non-thermal power-law indices yare near-identical (2.7 ±0.1), while 
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BBSO & SOHO/MDI &c HESSI,23-Jul-2002 00:28:45.000 




Fig. 2.3 RHESSI hard X-ray images of SOL20()2-()7-23T0():35 (X4.8) (white thin contours) taken at 
00:28 UT, overlaid on Michelson Doppler Imager (MDI) neutral-line magnetograms (white thick contours) 
and (negative) Ha taken at 00:28:45 UT. One e xtended (coronal) and three compact (footpoint) hard X-ray 
sources are evident. From lZharkova et all i2005h . 



the spectral index of the non-thermal emission for the coronal source is substantially steeper 
(y=5.6±0.1). 



2.3.2 Hard X-ray and gamma-ray sources 

Gamm a-ray ima ges were first obtained by RHESSI for the flare SOL2002-07-23T00:35 
(X4.8) ( IHurford et al. 2003) . They revealed the striking result that the hard X-ray sources are 
separated by a considerable distance (>20") from the 2.223 MeV neutron-capture gamma- 
ray source. This presumably reflects differences in the acceleration and/or transport pro- 
cesses for ions and electrons, respectively (Vilmer et al. 2011). A similar, but sma ller, sepa- 
ration was also detected for the "Halloween" flare SOL2003-10-28T1 1:10 (X17.2) faurford et all 
l2006h . 



12 



Zhai'kova et al. 




'0 100 10 100 10 100 

Energy (keV) Energy (keV) Energy (keV) 



Fig. 2.4 Comparison of the spectra of the footpoint and coronal sources of SOL2005-07-13T14:49 (M5.0). 
Left and middle panels: Footpoint spectra. A power-law has been fit between the dotted vertical Hnes. Right 
panel: Spectrum of the coronal source. A power-law plus an isothermal spectrum has been fit between the 
dotted vertical lines. The bold line is the power-law fit in the dashed interval and a Maxwellian fit for the coro- 
nal source at low-energies. The thin curve is the background (mostly instrumental) determined at a location 
outside the enhanced flare emission. From .Battaglia & Benz 1 2006) . 



2.4 Hard X-ray and radio emission 

The information obtained through study of radio emission from s olar flares complem ents 
that obtained from analysis of hard X-ray emission in several ways ( IWhiteetalJl201lh . Ra- 
diation at centimeter and millimeter wavelengths is gyrosynchrotron emission originating 
from mildly relativistic electrons with energies above about ~ 100 keV. Analysis of this 
emission hence permits an estimate of th e ambient magnetic field strength, reported to be 
~200 G in the corona jWhiteetal .11201 ll) . Coronal sources have been observed to emit in 
centimeter radio emission even before X-ray footpoints appear jAsai et al.l2006h . suggesting 
that the pre-flare heating may not be unrelated to acceleration in the main phase. Compari- 
son of radio and hard X-ray emissions bv iBenz et al.l ( |2003) also revealed, in some events, 
coherent radio emission after the hard X-ray emission had stopped (this phenomenon was, 
however, reported in only 5% of the flares studied). Gyrosynchrotron r adio emission in flares 
in general correla tes well with hard X-rays in time but not in space ( iKrucker et ani2008bl : 
IWhite et allioTlh : in addition, the electron spectra deduced from the radio emission are typ- 
ically harder than the electron spectra deduced from the hard X-ray emission. Flare emission 
at decimeter wavelengths concurrent with hard X-rays has also been found to originate from 
different locations: for example, Battaglia & Benz, (2009) report the origin of associated 
decimetric spike sources separated by some 100" from the coronal X-ray source. 

Emission at meter and decimeter wavelengths originates not from individu al particles, 
but fro m coherent radiation processes associated with various plasma waves. iBenz et al] 
ll2005h have studied the association of RHESSI hard X-rays with coherent radio emission 
in the meter and decimeter band in 201 flares with flare classes larger than C5. They found 
that coherent radio emission occurred before the onset of hard X-ray emission in 9% of the 
flares studied. Most of the radio emission was type III radiation at the plasma frequency and 
its first harmonic, emitted by a beam of electrons moving upward in the solar atmosphere 
into regions of progressively lower density. In a few cases, they also found a pulsating 
continuum, possibly caused by cyclotron maser emission. Both types of emission indicate 
the presence of electron acceleration in the pre-flare phase, which is often accompanied by 
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X-ray emission from a thermal source in the corona llKrucker et al.ll2008al : iFletcher et alJ 
l201lh . 

Emission at higher frequencies, such as decimetric narrow-band spikes, pulsations, and 
stationary Type IV events, correlate more frequently with th e hard X-ray flux and thus ap- 
pear to be more directly related to the acceleration process jArzner & Benzll2005l) . While 
there is a good association between coherent radio emission and hard X-rays, a strong cor- 
relation in the details of the time profile is less frequent, so that coherent emission is not 
a reliable proxy for the main flare energy release. This can be accounted for by invoking 
multiple reconnection sites connected by common field lines, along which accelerated par- 
ticles propagate and serve as a trigger for acceleration in other regions at some distance 
[fienz et aL^OOS.) . 

2.5 Magnetic field changes associated with the flare 

Sharp temporal increases of hard X-ray emission are often closely c orrelated in time with 
variations of the magnetic field measured on the photosphere (see, e.g.jKosovichev & Zharkoval 
I2OOII : ISudol & H arvev"2005l: IWang et aDl2005l : Izharkova et alJl20^ . For example, in the 
flare SOL2002-07-23T00:35 (X4.8), the magnetic flux change over the flare duration was 
about 1.2 X 10^^ Mx, most of which was concentrated in the area spanned by the magnetic 
inversion line. Such magnetic field changes are irreversible: the magnetic field reaches a 
new steady state and does not return to its pre-flare value (JCosovichev & Zharkova 200ll : 
ISudol & Harvevll2005l : Izharkova et alj|2005l) . Cross-correlation analysis with a time lag be- 
tween the temporal variations of magnetic field covering the neutral line and the hard X-ray 
light curve reveals a significant positive correlation with a time difference less than ~100 s. 
This observation strengthens the theoretical notion that irreversible changes in the magnetic 
field are responsible for the initiation and development of the flare phenomenon. 

3 Models of Particle Acceleration 

In this Section, we first review some basic concepts, such as acceleration by a direct electric 
field in a plasma environment (Section 13.1b and magnetic reconnection (Section 13.2b . in 
which magnetic energy is released into various forms, including the acceleration of charged 
particles, both electrons and ions. Recognizing the strong observational evidence for the 
formation of current sheets (localized regions in which a strong magnetic shear is present) 
in solar flares, our initial discussion relates to traditional two-dimensional (2-D) current 
sheets, with subsequent elaboration to the three-dimensional (3-D) case (Section[33]l. Two 
likely products of the fundamental reconnection process are magnetohydrodynamic shocks 
and stochastic MHD turbulence, and particle acceleration by these agents is discussed in 
SectionlSn 

3.1 Basic physics 

As remarked in Section[TJ the essential physics of particle acceleration is in principle rather 
straightforward - charged particles gain energy when subjected to an electric field in their 
rest frame. Such an electric field E may be a large-scale externally-imposed field, or a v x B 
field associated with particles crossing magnetic field lines, or a collective field associated 
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with the environment in which the particle finds itself (e.g., a coUisional Coulomb field or 
a field associated with a level of plasma wave energy). The richness of these various means 
of creating local electric fields is reflected in the richness of particle acceleration models in 

magnetized plasmas. 

One of the most basic concepts is that of the Dreicer field ( IChenlll974) . Consider an 
electron subject to an externally-imposed large-scale electric field E, plus the frictional force 
due to collisions with stationary ambient particles. The equation of motion for such a particle 
is (in one dimension) 

^=e£/m-/,v, (3.1) 

where m and e are the electron mass and (absolute value of) charge, fc is the collisio n 
frequency, and v is the velocity. Since for Coulomb collisions ~ v^"' (e.g., IChenll 19741) . 
we may write this as 

J4=e(E-Er>M-), (3.2) 



dt \ L V , 

where Vth is the electron thermal speed and we have defined the Dreicer field by 

E^^isM^, (3.3) 
e 

numerically, Ed ^ 10'^ n(cm"^)/7'(K) V cm"'. 

If £ < Ej), then only high-velocity particles (v > v„ii = Vth[£D/£^]'^^) see a net pos- 
itive force and gain velocity. This increased velocity reduces the drag force, leading to 
an even greater net acceleration and eventual runaway acceleration. In contrast, particles 
with V < Vcrit suffer a net retardation force and are not accelerated; they remain part of the 
coUisionally-redistributed Maxwellian core. For a Maxwellian distribution of velocities, the 
number density of accelerated particles /^^ /(v) dv (where /(v) is the phase-space distribu- 
tion function) is a (rapidly) increasing function of E. If E > Eo, all particles with velocity 
in excess of the thermal speed Vth are accelerated, and the efficient coUisional re-population 
of the accelerated electrons thus leads effectively to acceleration of all the electrons in the 
distribution. 



3.2 Magnetic reconnection models associated with flares 

It is generally accepted that the energy release in solar flares occurs through a reconstruction 
of the magnetic field, caused by the change of connectivity of magnetic field lines during a 
magnetic reconnection. The electric field associated with this changing magnetic field and/or 
with the associated driven currents leads to particle acceleration. 

Reconnection is a fundamental proces s defined by the magnetohydrodynamics of a mag- 
netized plas ma (see Priest & Forbesl2000h. The m agnetic diffusion equation (see, e.g., equa- 



tion (3.91) of lTandberg-Hanssen & Emslielll988h is 

^=Vx(vxB) + ^v2b, (3.4) 
at An 

where c is the speed of light, B is the magnetic field, v the fluid velocity and r\ the resistivity. 
The order-of-magnitude ratio of the terms on the right side of this equation defines the 
magnetic Reynolds number 

.„, = ^; (3.5) 
r\c^ 
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Fig. 3.1 Topology of the reconnecting magnetic field lines in the vicinity of an X-nuU point By = ±o;6^. 



it measures the ratio of the advective to diffusive contributions to the change in the magnetic 
field. For typical solar coronal values of rj and with £ w 10^ cm, R,„ w lO''*; for such high 
magnetic Reynolds numbers, the plasma is effectively "frozen-in" to the field and negligible 
change in the field topology, with its concomitant release of magnetic energy, can occur. A 
topological change in the magnetic field requires a breakdown in this ideal "frozen-in" flux 
condition; by Equation l3.5l this can occur either in small-scale regions (high values of V^B) 
and/or in regions in which the resistivity rj is anomalously enhanced. As a result, magnetic 
reconnection fundamentally occurs in narrow boundary layers called diffusion regions. 



3.2.1 Basic 2-D MHD theory of magnetic reconnection 

In the simple neutral sheet geometry originally proposed by Sweet and Parker jS weed 19691) . 
oppositely-directed magnetic field lines B in close proximity follow a plasma inflow V; ori- 
ented in the .^-direction perpendicular to the field lines - see Figure lSTl where the x-direction 
is vertical and the y-direction ho rizontal (note that there does exist evidence for such re- 
connection geometries in flar es - lAulanier eTal]|2000l ; iFletcher et al.ll200ll ;[ Sui & HolmanI 
l2003LlD"es Jardins et al.ll2009h . The increasing magnetic pressure in the localized region of 
field reversal is alleviated by reconnection in a small region of high V^B near the origin, 
which allows plasma outflo ws with velocity Vo toward the sides of the diffusion region. 
The original analysis bv lSweetl ( 1 19691) showed that the reconnection rate Vj is given by 

' (3.6) 



-1/2 

where Va is the Alfven velocity. Because of the dependence on the inflow veloc- 

ity, and hence on the rate of magnetic energy dissipation, such a neutral sheet scenario 
describes a "slow" magnetic reconnection regime. In an effort to increase the reconnec- 
tion rate. lPetschekl ( il964) proposed a model involving the formation of shocks propagating 
outward from the diffusion region; in this scenario, the reconnection process occurs at a 
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significantly faster rate, with V, w V^. Both these models were originally proposed with 
a classical value for the resistivity rj and both resulted in reconnection rates too slow to 
explain the rapid energy release observed in flares. However, later realization that the re- 
sistivity in the reconnection region may be driven to anomalously hig h values through the 
presence o f plasma instabilities in the current sheet (see, for example. iDahlburg et al.l[l992l : 
lBretl2009L and references therein) . These include the lower hybrid drift instability leading to 
the K elvin-Helmholtz instability (Landau & Lifshitz 1960; Chen et al. 1997; Lapenta e t al.l 
|2003|) . the kink inst abilitv (Pritchett et a l. 1996 ; Zhu & Winglee 1996; Daughton 1 999.), the 
sausa ge instability ( lBiichnejll999fr and magneto-acoustic waves ((Schekochihin & CowlevI 
l2007h . This realization revived hopes that reconnection mechanisms could indeed be a vi- 
able candidate for the energy release in flares (Svrovatskii 1981; Priest 1982; Somov 2003). 
Questions remaining about the stability of current sheet geometries such as those envisioned 
by Sweet and P arker have been extensively investigated using a full kinetic approach (see 
IPaughton et al.l l2009. and references therein), and are discussed further in Section l4~8] 

Another simple 2-D diffusion region geometry involves magnetic field lines of the form 
By = a}B\. Such a form is associated with two types of neutral point; an 0-type, ellip- 
tic topology {a} < 0; see Section l4!9l l. and an X-type, hyperbolic topology (a^ > 0; see 
Section [331 Figure [m . 



3.2.2 Effect of the Hall current on reconnection rates in 2-D models 

In the conductive plasma within the diffusive region where reconnection occurs, the current 
sheet tends to diffuse outwards at a slow rate with a timesca le of Tdif = a^S? wher e a is 
the current sheet half-width and fjj is the magnetic diffusivity jPriest & Forbesll2000l) . This 
diffusion process ohmically converts magnetic energy into heat. As discussed in Section [33] 
this timescale is too slow to account for the fast reconnection rates observed in the impulsive 
phase of flares unless some other processes, such as resistive instabilities or Hall current 
effects, are invoked to enhance the reconnection process. The effects of the Hall current in 
MHD models of the reconnection process are now discussed. Its effects in particle-in-cell 
(PIC) models will be discussed in Section l4~8l while the effect of resistive instabilities will 
be discussed in Section l4~9l 

The basic Sweet-Parker 2-D reconnection model discussed above can produce more 
realistic reconnection rates if, instead of a resistive MHD model, one uses a two- fluid MHD 
model, which includes in the equation of motion a Hall current ('Sonnerupl ll979h . Use of a 
two-fluid model means that the velocity vector v in Equation l3.4l now represents the velocity 
V of the plasma as a whole, with an additional vector V// representing the relative velocity 
between the electrons and the ions. The Hall current results from the fact that, inside the 
current sheet, the electrons and ions are magnetized to significantly different degrees - the 
spatial scales for effective magnetization is c/Wpe for electrons and c/Wpi for protons. (Here 
Wpe and (Opi are the respective plasma frequencies, which are in the ratio of the square root 
of the proton-to-electron mass ratio ^Jrn^Jme 43.) 

At length scales between c/Wpe and c/Wpi, the motions of the electrons and ions are de- 
coupled, and the current generat ed from the relative motion of electrons and ions can lead to 
the formation of whistler waves jMandt et alJl994d ; lMr& Bhattacharieell996l ; lDrake et alj 
ll997l ; lBiskai^ll997l) . Inside such a region, the reconnection rate becomes insensitive to 
electron inertia, allowing the reconnection rates to be stron gly increas ed to levels large 
enoug h to account for those measured in solar flares ( iBirn et ani200ll ; iHuba & Rudakovl 



l2004h . These results indicate a compelling need to consider a full kinetic approach for mod- 
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Fig. 3.2 The magnetic field geometry formed by four point-like sources of different magnetic polarity. The 
dot-dashed line labelled 'NL' is the neutral line. The separatrix surfaces are shown as thin dome-like lines 
conn ecting the photospheric sources. The intersection of these surfaces is the separator, shown as a solid line. 
From lTitov & Homid I2003) . 



eling magnetic reconnection processes in the diffusion region (for details, see Bim et al. 
l200ll . and references tiierein); this will be discussed in more detail in Sections [4 . 8 1 and l4!9l 



3.2.3 Further improvement of reconnection models 

It must be recognized that magnetic reconnection is inherently a 3-D process, and that the 
addition of the third dimension inherently introduces fundamentally new physics. (As an ele- 
mentary example, 2-D models may predict unrealistic energy gains since the induced electric 
field E = (V/c) X B is infinite in extent in the invariant direction). In three dimensions, the 
magnetic field has a considerably more complicated structure (Figure [Til l, involving sep- 
aratrices and separators (surfaces separating different domains of magneti c connectivity , 
and lines of intersect i on of separatrix surfa ces, respective lv; see Priest & Demoulinlll995l : 
iDemoulin et al.ll99dlbemoulin et al.ll99d) . As shown bv lPriest & DemoulirJ ( Il995h . mag- 



netic null points are no longer required for reconnection to take place, and, even when they 
do exist, magnetic nulls have a more c omplicated structu re, typically involving a "fan" struc- 
ture oriented around an axial "spine' ' jPriest etal.ll 19971) . 

In summary, reconnection can proceed in a variety of different magnetic field geometries 
and topologies (Priest & Forbe s 2002). Furthermore, depending on the magnetic field con- 
figuration, recormection can be accompanied by a variety of other physical processes such 
as generation of plasma waves, turbulence and shocks, each occurring on different temporal 
and spatial scales. These additional processes not only directly accelerate particles (e.g., in 
regions of enhanced magnetic field, which act as moving "mirrors"), but can also feed back 
on the magnetic field topology and hence the ongoing reconnection rate. 
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Plasma inflow/outflow x 




Fig. 3.3 Top: current-sheet location and bottom magnetic-field topology in a single helmet-type reconnecting 
non-neutral current sheet. The large arrows depict plasma inflows from the top and the bottom and outflows to 
the side. The x- and z-components of the magnetic field lie in the plane of the Figure; a guiding magnetic-field 
component By and a drift electric field are perpendicular to the plane. 



3.3 Particle acceleration in a reconnecting current sheet 

Two of the most popular ma gnetic field topologies associated with solar flares involve the 
formation of a current sheet ( [Gorbachev & Somovll 19891 : ISui & Holmai]||2003l) , located ei- 
ther at the top of single helmet-like loop (Figure [33] | or at the intersection of two interacting 
loops (Figure [34t . Other more complicated topologies are essentially combinations of these 
two basic ones, and they are discussed in Section |4] In this Section we therefore briefly 
overview the basic physics of particle acceleration in reconnecting current sheets. 
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Fig. 3.4 Top: cuiTent sheet location, and bottom: magnetic-field topology in a reconnecting non-neutral cur- 
rent sheet formed by two interacting loops. The large arrows depict the plasma inflows. As for Figure [331 the 
X- and ^-components of the magnetic field lie in the plane of the Figure; a guiding magnetic field component 
By and a drift electric field are perpendicular to the plane. 



Let us, then, consider the magnetic field topology as presented in the reconnecting cur- 
rent sheet (RCS) geometries illustrated in bottom plots of Figures l3.3l andl 3.4l The tangential 
magnetic field component B-^ {x) is approximated by the linear function 



B,{x)=Bo (-^), 



(3.7) 



where d is the RCS thickness jMartens & YoundliggOlJLitvinenko & Somovlll993bl) . The 
width a of an RCS in the z-direction approaches 10^ - 10^ m i Somov[l992 ) and magnitudes 
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of the other two magnetic field components Bx and By are assumed to be constant: 

|B.,|=S,<,, |Sv|=S,,o. (3.8) 

Particles, effectively frozen into the field lines, are advected toward the RCS plane (the 
X-axis) by the flow field V,, where they are accelerated by the (Vj/c) x B e lectric field 



x-axis) by trie tlow tield Vj, wnere tney are accelerated by trie (Vj/c) x j> e lectric tielc 
teulanov & Svrovatskii|[l97^ : iLitvinenko & Somovlll993bl : lLrtvinenkdll996bl) . As the re 
connection process continues, ambient particles drift into the reconnecting current sheet; the 
rate of acceleration of particles is limited b y the size of the thin diffusion region. 

T he induced electric field is given by i Bulanov & Svrovatskiill97^ : lMartens & Yound 



1 1990l : ILitvinenko & Somovlll993b| - |zharkova & Gordovskvvll2004h 



Ey^Y^-f'-^; (3.9) 
c An ax 

the first and second terms on the ri ght hand side o f Equat i on|3.9l dominate outside and inside 
the diffusion region, respectivelv ('S omovll2000l) . IPriestI ( Il982) has shown that the plasma 
inflow velocity Vi is smaller than the local Alfven speed by about two orders of magnitude, 
i.e., Vi w O.OIVa- Hence the electric field is 

252 

0.01^=1^, (3.10) 
c ^JAnnnip 

where nip is the proton mass. Let us assume the following physical conditions in the current 
sheet: temperature T = 10* K, density n = 10'° cm"^ and B-q w 100 G. Then Va~2x 
10* m s^' and the thermal velocities are about 10^ m s^' for protons and 10^ m s^' for 
electrons. Taking the initial inflow velocity V; to be the proton thermal speed, we obtain 
Z?v ~ 2 V cm^ ' , which is rou ghly four orders of magnitude higher than t he corresponding 
Dreicer field -lO""* V cm"> ^Martens & Youndll990l : lLitviiienkdll996bl . see Section [37B. 
This electric field is directed along the longitudinal magnetic field component B,,, allowing 
efficient particle acceleration in the y-direction. 

In these simulations the drifted, or reconnection, electric field is assumed to exist through- 
out the entire diffusion region, which in the MHD approach extends over som e 10^ m to 
10 m teulanov & Svrovatskii|[T97^ : ISomovlll992l : [ Oreshina & Somov ''2000a'). This as- 



sumption has been challenged by full kinetic PIC simulations ( Furth et al.ll963i : , Drake & Led 
Il977h . which show the reconnection field to extend over a relatively small domain corre- 
sponding to the size of the electron diffusion region. Given the current limitations of PIC 
simulations (see Section 14.10b . there is at present no definitive indication of which of the 
approaches (MHD or PIC) is more valid. In the remainder of this Section we adopt the 
value provided by the MHD approach, recognizing that appropriate caveats must therefore 
be made. 

The magnitudes of the magnetic field components B^ and By are chosen such that elec- 
trons remain magnetized in the vicinity of the midplane while protons are unmagnetized 
during the acceleration phase . In this case, the energy of the unmagnetized protons upon 
ejection can be estimated (see lLitvinenkd[l996bl) as: 

e = 2mac'(^\ , (3.11) 



vO 



where nia is the particle mass. Howev er, since the electr ons are magnetized, they follow 
magnetic field lines and gain an energy ( lLitvinenkdri996bh : 



eEyo — 
Bxi) 



(3.12) 
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where d is the current sheet thickness along the x-direction (see the bottom plot in Figure 
I3.3I >. For RCSs with thicknesses from 1 m to 100 m, corresponding to proton Larmor radii 
in magnetic fields from 10^^ T to 2 x 10^^ T, electrons gain energies up to ~100 keV, 
comparable with those of bremsstrahlung-emitting electrons. 

Although such a model does meet the challenge of (1) a sufficiently high reconnec- 
tion rate and (2) sufficient energy gain by the accelerated electrons, the number of accel- 
erated particles is fundamentally limited by the small width 2a of the current sheet. These 
deficiencies are addressed in the further development of this basic model to include three- 
dimensional considerations (Section [4.41 ). the use of multiple current sheets formed by the 
filamentation of the reconnection region (Section 14.6b , and a full kinetic PIC approach to 
model portions of a 3-D current sheet (Sections 14. 81 and |4~9l l . 

3.4 Particle acceleration by shocks and turbulence 

In this Section, we review the basic physics involved in acceleration models that do not rely 
on a large-scale coherent electric field. 

3.4.] Shock acceleration 

Ever since the demonstration by several authors telandford & Ostrike j [ToTi : lBdilll978l ; 
ICargill 1991 ; Giacalone & Joki pii 1996) that a very simple model of shock wave accelera- 
tion leads to a power-law spectrum of the accelerated particles (in good agreement with ob- 
servations of galactic cosmic rays), shock acceleration has been frequently invoked in both 
space and astrophysical plasma contexts. However, this simple model, though very elegant, 
has a significant number of shortcomings. Inclusion of losses (Coulomb at low energies and 
synchrotr on at hi gh energies) or the influence of accel erated particles on the shock structure 
(see, e.g.. lAmato & Blasi 2005: Ellison & et alj|2005h can cause significant deviations from 
a simple power law. 

Moreover, shocks are unable to accelerate low energy background particles efficiently; 
shock acceleration requires the injection of (fairly high energy) "seed" particles. Also, shock 
acceleration requires some scattering agents (most likely plasma waves or turbulence) to 
cause repeated passage of the particles through the shock front. In the case of quasi-parallel 
shocks, in which the magnetic field is nearly parallel to the normal of the shock front, the 
rate of energy gain is governed by the pitch-angle diffusion rate. 

As pointed out by Wu (1983) and|Lerov & Mangeney (1984), in the Earth's magne- 
tosphere, many shocks can be quasi-perpendicular (i.e., the vector defining the normal to 
the shock is almost perpendicular to the shock velocity vector). In such a case the par- 
ticles drift primarily along the surface of the shock and get accelerated more efficiently. 
However, acceleration in quasi-perpendicular shocks suffers from the fundamental limita- 
tion that injection of seed particles with very large velocities v > u^h^ (where is the shock 
speed and £, ^ I is the ra tio of parallel to perpendicular diffusion coefficients) is required 
llAmano & Hoshindl2010h . Moreover, although there may b e some indi rect morphological 
evidence for the existence of such shocks in flares (see, e.g.. lSui & Holman 2003), there is 
little direct evidence for the occurrence of these kinds of shocks near the top of flaring loop, 
where the acceleration seems to be taking place during the impulsive phase. Furthermore, 
such shocks typically do not appear in 3-D MHD simulations of the reconnection process. 

In summary, many of the features that make acceleration of cosmic rays by shocks at- 
tractive are not directly applicable to the solar flare environment. As a result, this mecha- 
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nism has, with some justification, received only limited attention in the literature. However, 
in Section 14.3.21 we do discuss the possibility of acceleration at fast magnetosonic shocks 
resulting from reconnection. 

3.4.2 Stochastic Acceleration 

A main contender for the particle acceleration mechanism in flares is stochastic acceleration 
by the magnetic field component of low-frequericy ma gneto-acoustic waves jMiller et alj 
ll997l : ISchlickeiser & Mille3ll998l : |Peb^osian et al.ll200d) . 

In the originally proposed Fermi acceleration process ('F ermilll949l : lDavis| [ 195a), par- 
ticles of velocity v scatter coherently off agents moving toward each other with a veloc- 
ity ±M. Since the velocity gain in a single collision is 2u, the rate of velocity increase is 
Av/ At w 2m/(L/v), where L is the separation of the scattering centers. Solving this, for the 
case in which a scatter occurs on a timescale significantly less than that required for L to 
decrease, shows that the energy grows exponentially with time: E cx g^u'//. p^^ particles 
scattering off randomly moving agents, both energy gains and losses are possible. However, 
the increased likelihood of a head-on (energy-increasing) collision leads to a general (slower 
than exponential) increase in energy, at a rate proportional to [u/v]^ D^j^, where D^^i is the 
angular diffusion coefficient. This process is known as second-order Fermi acceleration. 

In general, in a magnetized plasma characterized by a gyrofrequency Q , a particle of 
velocity v will resonate with a wave of frequency CO and wavenumber k when 

a = k\-m, (3.13) 

where i is an integer. Particles with velocities near the Cherenkov resonan ce (v w cq/A: || ) 
resonate with i = waves. This process is known as transit-time damping iHollweg||l974 . 



and references therein). It acts to diffuse the particle distribution /(p) in momentum space, 
described by the Fokker-Planck equation 

The diffusion coefficient Djj and the secular coefficient (; = 1 , 2 correspond to energy and 
angular terms, respectively) contain the essential physics: the action of both the accelerating 
waves and decelerating effects (e.g.. Coulomb collisions); the quantity 2(p) corresponds to 
the injection of particle energy and the term 5(p)/(p) to escape. 

For solar fl ares, the most likel y sources for the energization t erm 2(p) are plasma wave 
turbulence ([Hamilton & Petrosiarj|l992i : IPrvadko & Petrosianlll997l) and cascading MHD 
turbulence i Miller et alT 1996h . In the first scenario, a presupposed level of plasma wave 



turbulence in the ambient plasma accelerates am bient particles of energy E stochastically 
with the rate Dee/E^- It has been shown (e.g.. iHamilton & Petrosianlll992l : iMiller et alj 
1 19961) that for flare conditions, plasma-wave turbulence can accelerate the background par- 
ticles to the required high energies within the very short timescale derived from hard X- 
ray observations. More importantly, at low energies, and especially in strongly magnetized 
plasmas, the acceleration rate Dee/E^ of the ambient electrons may significantly exceed 
the scatt ering rate Dug for particula r pitch angles and directions of propagation of this tur- 
bulence dPrvadko & Petrositujll997l . Il998h . leading to an efficient acceleration. However, 
high-frequency waves near the plasma frequency are problematic as drivers for stochastic 
acceleration in flares, because they would couple into decim eter radio waves a nd would be 
present at high levels in every flare, contrary to observations ( iBenz et al.ll2005h . 
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In the second scenario, involving cascading MHD turbulence, moving magnetic com- 
pressions associated with fast-mode MHD waves propagating at an angle to the magnetic 
field constitute the scattering centers. The wavelength of the induced MHD turbulence is 
assumed to cascade through an inertial range from large to very small, stopping at the point 
where wave dissipation becomes significant due to transit-time damping. Resonance with 
one w ave will result in a n energy change leading to resonance with a neighboring wave, and 
so on dMiller et 31.11199 6*). Electrons can thus be accelerated all the way from thermal to rel- 
ativistic energies through a series of overlapping resonances with low-amplitude fast-mode 
waves in a continuum broad-band spectrum. 

In order to be efficiently accelerated, the test particle must have an initial speed that 
is greater than the speed of the magnetic compression, i.e., of order the Alfven speed. For 
the pre-flare solar corona, electrons have speeds comparable to, or greater than, the Alfven 
speed, and so are efficiently accelerated immediately. On the other hand, in a low-jS plasma, 
the protons necessarily have an initial speed that is less than the Alfven speed and so are not 
efficiently accelerated: the scattering centers converge before the test particle has suffered a 
sufficient number of energy-enhancing collisions. Efficient acceleration of ions therefore re- 
quires pre-acceleration to Alfven speeds. In the model of Miller & Roberts ( 1995J, this pre- 
acceleration is accomplished through (parallel-propagating) Alfven waves, and it is essential 
to realize that this process takes some time to accomplish. During this pre-acceleration pe- 
riod, most of the energy is deposited in accelerated electrons; only after the protons are 
pre-accelerated to the Alfven speed are they effective ly accelerated - they the n gain the bulk 
of the released e nergy by virtue of their greater mass dMiller &Robertsll 19951) . As explained 
bv lEmslie et al.l (2004). this may be the explanation for the different locations of h ard X-ray 
and gamma-ray sources in flares (see Section |23] and iHurford et al.1l200ll2006h : efficient 
proton acceleration requires relatively long trapping times and so a larger magnetic structure, 
whereas electrons can be efficiently accelerated even in relatively small structures. 

In practice, since there is a lit tle differenc e mathematically between the plasma wave 
turbulence and shock mechanisms dJ onesll 19941) . acceleration by both turbulence and shocks 
can be effectively combined. It should be also be noted that stochastic acceleration generally 
predicts an accelerated spectrum that is not a strict power-law, but rather steepens with 
mcreasing energy fe.g.. lMilleret 'ZI [l99^ . However, a recent development in transit-time 
damping theory involves adding both the loss term ^(p) in Eauation l3. 14l that allows electron 
escape, and a source term 2(p) due to the return current replenishing the electrons in the 
acceleration region. These additions c ause the solution to be come stationary within less than 
one second under coronal conditions dGrigis & BenzlbOOd) : furthermore, in the hard X-ray 
spectral region of observed non-thermal energies, the photon distribution is close to a power 
law, in agreement with observation. The predicted spectral index y depends on the wave 
energy density and on the interaction time between the electrons and the waves. Additional 
trapping (e.g., by an electric potential) tends to harden the electron spectrum. 

The main limitation of stochastic acceleration models at this juncture is the need to 
presuppose an ad hoc injection of the necessary waves (plasma or MHD). Furthermore, in 
the absence of pitch-angle scattering, the stochastic acceleration process generally leads to 
a decrease in particle pitch angles. This, in turn, reduces the efficiency of the acceleration 
process, since now only waves with very high parallel phase speeds, or those with pitch 
angles near 90°, can resonate with the particles being accelerated. Therefore, in order for 
significant stochastic acceleration to occur, the effective particle scattering rate (which is 
proportional to the collision frequency between particles and waves) must be greater than 
the frequency of Coulomb collisions between the charged particles. The latter condition 
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imposes rather strict limitations on the growth rates of plasma waves or MHD turbulence, 
as appropriate. 



4 Recent Theoretical Developments 

In this Section, we highlight some recent promising developments in our understanding of 
the physics of particle acceleration in the magnetized plasma environment associated with 
solar flares. 

In Section l4n we turn our attention to some recent advances in stochastic acceleration 
models, namely the inclusion of (1) small-scale electric fields and (2) feedback of the accel- 
erated particles on the turbulence level in the acceleration region. Some ideas on possible 
particle acceleration by large-scale Alfven waves in a large-scale closed magnetic configu- 
ration are discussed in Section l4~2l 

In Section |43] the evolution of the ubiquitous "helmet"-type magnetic field configura- 
tion observed in flaring regions is presented. We point out two likely acceleration regions in 
such a geometry: the collapse of previously-stretched field lines onto a closed magnetic loop 
(which leads to acceleration by both first-order Fermi and betatron processes, associated 
with the shrinking and compression of field line bundles, respectively), and the formation 
of "termination" shocks associated with high-speed jets, which accelerate particles through 
shock-drift acceleration. 

In Section l44l we consider particle trajectories in a realistic 3-D current sheet geometry, 
considering, in particular, the asymmetric acceleration of electrons and ions and the result- 
ing particle energy spectra, which turn out to be power laws. In Section |431 we discuss 
particle acceleration in 3-D magnetic field configurations, in which magnetic nulls have a 
fundamentally different topological structure than they do in a simple 2-D reconnection ge- 
ometry. Both "fan" and "spine" reconnection modes are considered, and the resulting (rather 
hard power-law) spectra of the accelerated particles are presented. 

In Section l4~6l we turn our attention to acceleration in multiple acceleration sites, con- 
sidering cellular automaton (CA) models, MHD turbulent environments, and stochastic ac- 
celeration in a distribution of elementary current sheet sites. 

In Section l477l we turn our attention to a critical aspect of all acceleration models - the 
need to recognize that, since a substantial fraction of the magnetic energy released in flares 
appears in the form of accelerated particles, these particles simply cannot be considered 
as an ensemble of test particles interacting with a prescribed electrodynamic environment. 
Rather, the electric and magnetic fields generated by the accelerated particles themselves 
are an essential element of this environment; the self-consistent inclusion of these fields 
represent a paramount consideration in acceleration models, without which unphysical and 
even paradoxical conclusions can result. We review the limitations of existing acceleration 
models in this context. In Section l4~8l we discuss recent attempts to include this electrody- 
namic feedback through the use of PIC models, for which the self-consistent electric and 
magnetic fields associated with the accelerated particles do form an essential element of the 
modeling. Finally, in Section l4~9l we briefly discuss a mechanism for particle acceleration 
that involves multiple small-scale closed magnetic "islands" formed by the tearing-mode 
instability in a reconnecting current sheet. 
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4. 1 Stochastic acceleration 

As discussed in Section IZTl tlie low value of the low-energy cutoff of the electron energy 
distri bution suggested by RHESSI observations jKontar et al.ll200j : ISaint-Hilaire & Bend 
l2005h has significant implications for the number of accelerated flare electrons. Indeed, 
for sufficientl y low cutoff energies, the number of electrons required in a coUisional thick- 
target model te rownlligTll) becomes so great that models in which the accelerated particles 
form a subset of the ambient distribution become untenable. The only viable mechanism of 
particle acceleration is then acceleration of the entire plasma within the acceleration region 
through a stochastic process that results in little or no net electrical current, to avoid the 
issues presented in Section lZTl 

As discussed in Section [T4l high-frequency waves near the plasma frequency are prob- 
lematic as drivers for stochastic acceleration, because they would couple into decimeter 
radio waves and would be present at high levels in every flare, contrary to observations 
jBenz et al. ![2003). On the other hand, acceleration by transit-time damping do not suffer 
from this drawback - the driving waves have a frequency that is far below the plasma fre- 
quency. 

4.1.1 Improvements to the MHD stochastic acceleration model 

The basic transit-time damping stochastic acceleration mechanism has been improve d by 

addin g low-frequency fluctuating electric fields parallel to the magnetic field ( Grigis & Bena 

I2OO 4|,|2005). Such fields may originate from low-frequency and high-amplitude turbulence, 

such as kinetic Alfven waves, as discussed in Section lSTI fsee also lMiller et al.ll996l : |Prvadko &. PetrosianI 

Il997h . They can accelerate and deceler ate electrons (and protons), leading to a net diffusion 

in energy space jArzner & Vlahosl2004h . As for other models for stochastic acceleration, the 

authors assume an ad hoc distribution of turbulence. In this model, the non-thermal electron 

distributions in coronal sources "grow" out of the thermal population: the coronal source is 

initially purely thermal, then a soft non-thermal population develops, getting harder at the 

hard X-ray peak and softening toward the end of the emission. 

As discussed in Section 12. II such a " soft-hard -soft" behavior has b een observed in 
a coronal source ( iBattaglia & Benj|2007t) . Indeed, iGrigis & Be"n3 ( |2004) report RHESSI 
observations suggesting the existence of a point in the spectrum at which the hard X-ray 
flux does not change with time during an event. Above this pivot point the spectral varia- 
tion is consistent wi th the characteristic "soft - hard-s oft" behavior of the photon spectrum. 
Note, however, that IZharkova & Gordovskwl ( |2006|) have pointed out that such a "pivot- 
point" behavior in the variation of the hard X-ray spectrum with time could also be associ- 
ated with an electron beam accelerated by the super-Dricer electric field in a current sheet 
and precipitating into a l oop leg with the self-induced (return-current) electric field (e.g., 
iKnight & Sturrockll 1 977l : lEmslidl 198 ih associated with an electron beam accelerated by a 
directed electric field. 

4.1.2 Stochastic acceleration with particle feedback 

A mor e self-consistent approach to stochastic acceleration was recently reported bv lBvkov & FleishmanI 
(l2009h . who simulated electron acceleration by MHD waves, but included the effect on the 
MHD turbulence (spectrum and intensity) caused by the injection of high energy particle 
beams into the acceleration region. The initial MHD turbulence is prescribed, and is as- 
sumed to be produced in the form of transverse motions with a Gaussian spectrum on the 
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Fig. 4.1 The reconfiguring coronal field launclies a torsional Alfven wave pulse thr ough the corona and into 
the chromosphere, as well as a fast-mode wave pulse. The Alfven wave, which propagates in the inertial 
regime, can lead to electron acceleration in the corona. Any Alfven wave energy entering the chromosphere 
can also lead to electron acceleration there. The wa ve will be partially reflecte d from the steep gradients in 
the chromosphere and then re-enter the corona (from lFletcher & Hudsonl2008l) . 



scale 27z/ko. Through mode coupling, this turbulence produces a corresponding level of 
longitudinal turbulence in a system of a finite size R (where > !)• The model includes 
only large scale energy-conserving motions, with the kinetics of particles on smaller scale 
being determined by turbu lent advection, valid only for electron energies up to ~1 MeV 
tevkov & Fleishmarjl2009h : thus, the model can account for the feedback of the accelerated 
particles on the acceleration mechanism on large scales only. The phase-space diffusion co- 
efficients Dij in Equation [3TT4] are expressed in terms of the spec tral functions that describ e 
correlations between large scale turbulent motions (for details, see lBykov & To ptvgin'l993h. 

The model assumes a continuous injection of monoenergetic particles (both electrons 
and protons) into the acceleration region. The initial phase of the acceleration is charac- 
terized by a linear growth regime which results in effective particle acceleration by the 
longitudinal large scale turbulent motions and thus leads to a spectral hardening. However, 
because the accelerated particles eventually accumulate a considerable fraction of the tur- 
bulent energy, the efficiency of the acceleration decreases; this predominantly affects the 
higher-energy particles, leading to a spectral softening. 



4.2 Particle acceleration by large-scale Alfven waves 

In this model, reconnection is still at the heart of the magnetic energy release process, in that 
it allows the pre-flare field to reconfigure, releasing the energy needed to power the flare. 
However, motivated by in situ observations and theories for e nergy transport and ele ctron 
acceleration in the Earth's ma gnetosphere during sub-storms jStasiewicz et al.ll2000l) . it is 
suggested dPletcher & Hudson! 12008) that the energy is transported through the corona not 
by particle beams but by Alfven-wave pulses, launched in the primary reconnection site and 
carrying a large Poynting flux towards the ch romosphere. The presence of such pulses has 
recently been verified in MHD simulations bv lBirnetal.l ( l2009h . They may play a direct role 
in accelerating the electrons either in the field-aligned electric fields that they generate in the 
corona, or through the magnetic energy that they transport to the chromosphere, leading to 
acceleration of electrons locally in the chromosphere by mechanisms as yet unspecified. A 
cartoon realization of this model is shown in Figure |4~T] 

To explain the rapid time variability observed in the flare time profiles, and also the 
simultaneity of hard X-ray footpoint sources within around 0. 1 s, the wave pulses must move 
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quicldy - on the order of 0.1c, which is significantly larger than the values usually quoted 
for the coronal Alfven speed. However, such speeds are consistent with the properties of the 
core of an active region before a flare, where the density ma y be 10^ cm~-^ or less , and the 
magnetic field on the order of a few hundred Gauss (e.g.. iBrosius & White! 1 200^ . Linear 
wa ve pulses traveling a t such s peeds do not damp significantly en route through the corona. 

iFletcher & HudsonI bOOSh calculate that for a coronal field of 500 G and an ambient 
density of 10^ cm^^, a magnetic field perturbation of 50 G is more than enough to de- 
liver sufficient Poynting flux to the chromosphere to produce the flare energies observed. 
Of course, such a perturbation, on the order of a few to 10%, does stretch the linear ap- 
proximation for Alfven waves. An ideal, incompressible, non-linear Alfven mode traveling 
in a com pletely uniform plasma do es not steepen into a shock, or dissipate due to Landau 
damping jSames & Hollwedl974) . but the presence of density fluctuations, including those 
generated by the wave itself, can lead it to decay via a parametric instability into a for ward- 
propa gating ion-acoustic wave and a backward-propagating Alfven wave (e.g. Golds teiril 

Analytical work bv lJavanti & HoUwed ( fl993h suggests that, in the low-/3 coronal regime, 
the maximum growth rate of this instability for a wave with an amplitude of a few percent 
is a few times the wave frequency, so that this decay occurs for all but long-wavelength per- 
turbations. However, in the low-;S case, the Alfven speed exceeds the sound speed, so the 
front of the Alfven wave train (and presumably also the front of a simple Alfven wave pulse) 
propagates ahead of the ion-acoustic wave produced, and is, therefore, not significantly af- 
fected by the instability. This has been demonstrated by Turkrnani & Torkelsson ( 200^) in 
numerical simulations for roughly the parameter regime of interest. The process has not yet 
been fully explored in an environment where coUisionless damping is important, though this 
may be the case for coronal depths in flares. 

Because the Alfven speed in such an environment is so high (higher than the electron 
thermal speed, such that an average electron lags behind the wave), the details of the wave 
propagation require two-fluid or eve n full kinetic s imulations. Such simulations performed 
in the case of the magnetosphere ( Kletz in§j|l994l) indicate that a parallel electric field is 
generated by the wave pulse, in which ambient electrons can - if at the appropriate "injec- 
tion" energy - be accelerated to twice the Alfven speed. Similar considerations for solar 
coronal conditions indicated that modest non-thermal energies of a few tens of keV could be 
achieved in this manner. .McClements & Fletcher (2009) have also recently confirmed this 
for solar plasma in the test-particle limit. This process is found to damp the wave energy, 
though it remains to be seen by how much. 

On arriving at the steep gradient in Alfven speed represented by the chromosphere, it is 
to be expected that this wave pulse would be at least partially reflected jEmslie & Sturroc3 
fl982i) . However, any component of the wave energy entering the chromosphere could be 
locally damped - for example by ion-neutral coupling - leading to chromospheric heating. 
If one considers the wave pulse as imparting stress to the chromospheric magnetic field, 
then the multiple current sheets that t his might generate cou ld provide an environment for 
particle acceleration, as suggested bv iTurkmani et alJ ll2005h . The timescales are also such 
that a perpendicular MHD cascade could be established in the reflected wave field, provid- 
ing another opportunity for stochastic acceleration of the chromospheric electrons, if the 
acceleration timescale is less than the coUisional loss timescale as discussed in Section 

Significant theoretical work is yet required to demonstrate theoretically whether any of 
these options are viable, e.g., whether any chromospheric electron acceleration mechanism 
can reproduce the timing delays observed between high and low energy hard X-rays (e.g., 
lAschwanden et al]|l996h , which are normally interpreted as strong evidence for a coronal 
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Fig. 4.2 Sketch of the flare scenario. EP: erupting prominence; DR: diffusion region; SMSS: slow-mode 
standing shock; FMSS: fast-mode standing (also called termination) shock ; PFL: post flare loop. The grey 
shaded areas with the arrows show the outflow jets. From lMann et alj j2006l) . 



electron accelerator. It should also be noted that in the proposed model of energy transport, 
magnetic reconnection still plays a central role as the mechanism of energy release, and 
as the driver of the Alfven waves. This means that concurrent electron acceleration in the 
corona can also take place by "traditional" means (e.g., by electric fields in current sheets or 
by turbulent acceleration by various MHD wave modes). 



4.3 Electron acceleration in collapsing current sheets 

Figure I4. 2 | depicts a helmet-ty pe magnetic structure often observed in solar flares (see, for 



example, ISui & Holmanll2003h . in which a prominence is destabilized due to photospheric 



footpoint motions and rises upwards as an eruptive prominence. As the prominence rises, the 
underlying field lines are stretched, leading to the formation of a current sheet. We suppose 
that when the current exceeds a critical value, the resistivity is suddenly enhanced (say by 
plasma waves excited by various instabilities) and rapid magnetic reconnection begins in the 
diffusion region (DR), as shown earlier in Figure [BTI 

MHD solutions for the large-scale magnetic configuration plotted in Figure l4~2l allow us 
to investigate (see Section [l!2l i the reconnection process inside the diffusion region, charac- 
terized by plasma inflows with velocity Vb in the jc-direction (from the sides in Figure 14.21 
from the top and bottom in Figure 13. Il l and outflowing with velocities V\ in the z-direction 
(to the top and bottom of Figure l4~2l to the left or right sides of Figure [STTT l. The outflows 
are generally very fast and create a rather complicated structure in the outer region at the 
RCS edge. The plasma confined on the newly reconnected field lines will start moving with 
the velocity V\, both upwards (not shown in the cartoon), and downwards, where it will 
encounter a magnetic obstacle in the form of the underlying loop. If the distance to this 
obstacle is large, a collapsing magnetic trap (see the bottom plot in Figure I4.3l l is formed. 
Moreover, the downward outflow is moving into a region of generally higher density and 
hence lower Alfven speed. If a point is reached at which the outflow velocity V\ exceeds 
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Fig. 4.3 Schematic models of collapsing magnetic field lines occumng in MHD models of magnetic recon- 
nection with (left panel, dashed line) and without (right panel) a termination shock. From lSomov & Kosugil 
fr997i) . HTTCS denotes "high temperature turbulent current sheet." 



the local magneto-acoustic wave velocity, then a fa st MHD shock, termed a "termination 
shock," can appear ahead of this magnetic obstacle dSomov & KosugilfT997h . as shown by 
the dashed line in the upper panel of Figure |43] 

4.3.1 Acceleration in collapsing magnetic traps 

ISomov & Kosugil (fT997h . lKarlickv & Kosugil (l2004h and liogachev & Somovl ( l200l [20071) 

have carried out simulations of the acceleration of electrons in collapsing magnetic traps. 
Particles with an initial Maxwellian distribution are injected at random points of a frag- 
mented current sheet geometry obtained using a three-dimensional MHD simulation. 

As a result of its motion from the reconnection region towards the chromosphere, the 
length L of the magnetic trap decreases, leading to an increase in energy of the trapped 
particles via a fi rst-order Fe rmi process associated with conservation of the longitudinal 
invariant / ds ( lBenj|2002h . Also, as a result of the transverse contraction of the magnetic 
trap, particles are accelerated by the betatron mechanism. This transverse contraction can 
be described by the ratio b{t) = B{t)/Bo which changes from 1, when 5(0) = Bq to b„, = 
Bm/Bo, where S„, is the magnetic field in the magnetic mirrors. Thus the particle momenta 
P|| and p± change as follows: 

p\\=po\\r\ (4.1) 

where £ = (L/Lq), and 

P± = Po±b'^\ (4.2) 

Both the Fermi and betatron accelerations are accompanied by a change in particle pitch 
angle - a decrease (Fermi) or increase (betatron). When these two mechanisms act simulta- 
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neously, then the pitch angles of accelerated particle vary according to 

tana = — =^fc'/2 tan oo, (4.3) 

P\\ 

which can result in an increase or decrease of pitch angle, according to the relative impor- 
tance of the two effects. The particle kinetic energy increases in the magnetic trap until it 
falls into the loss cone. 

Simulations show that particles with an initially Maxwellian distribution gain energies 
via both acceleration mechanisms; they retain a quasi- thermal distribution but attain much 
higher energies. The change of total kinetic energy has a form of the burst of particles with 
energies up to a few MeV. However, the magnetic energy is absorbed very rapidly by the 
electrons, so that a test-particle approach very rapidly ceases to be valid. A fuller treatment 
of this scenario therefore requires full kinetic simulations similar to those discussed in Sec- 
tion Ej] 



4.3.2 Acceleration by a termination shock 



Due to the strong curvature of the magnetic field lines in the vicinity of the diffusion re- 
gion in a helmet-type reconnection (see Figure I4.2| |. the slowly inflowing plasma shoots 
away from the reconnection site as a hot jet. These oppositely-directed jets are embedded 
between a pair of slow magnetosonic shocks. If the speed of this jet is super- Alfvenic, a fast 
magnetosonic shock, also called a termination shock, is established, as shown previously in 
Figure 1431 The appeara nce of such shocks was predicted in the numerical simulations of 
iForbesI ( fl983) , [Shibata eTall j 19951) a nd Somov & Kosudl (Il997h . 



Radio measurements jAurass etal.,.2002. : .White et alfeoill) of SOL2003-10-28T11:10 



(X17.2) reveal the termination shock as a standing radio source with a (large) half-power 
source area of Aj = 2.9 x 10^" cm^. A strong enhancement of the electromagnetic emission 
in the hard X- and gamma-ray range up to 10 MeV is simultaneous ly observed with the 
appearance of the radio signatures of the shock jAurass & Maiinll 20041) . These observations 
indicate that the termination shock could well be the source of the highly energetic electrons 
needed for the generation of the hard X- and gamma-ray radiation. 

A fast magnetosonic shock formed during the MHD simulations of a reconnection pro- 
cess is accompanied by a compression of both the magnetic field and the density. Thus, it 
represents a moving magnetic mirror at which charged particles can be reflected and accel- 
erated through a process termed sh ock-drift accelera tion. This process can accelerate highly 
energetic electrons up to 10 MeV jMann et aDl2006h . 

Generally, shock-drift acceleration represents reflections at the termination shocks. Anal- 
ysis of such an acceleration is most straightforwardly performed in the de Hoffmann-Teller 
frame, the frame in which the electric field vanishes, so that the reflection takes place 
und er conservation of both kinetic energy and magnetic moment. From such an analy- 
sis jMann et"aD l2006'. 2009), one obtains the relationships between the electron velocities 
j3 = v/c parallel and perpendicular to the upstream magnetic field before (index /) and after 
(index r) the reflection: 

2jS,-ft||(l+j62) 



and 
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where jSj = sec d /c. Here Vj is the shock speed, c the velocity of Ught, and 6 the an- 
gle between the shock normal and the upstream magnetic field. In addition, the reflection 
conditions: 

ft.li < ft (4.6) 

and 

ft,,>-J^.(ft-A,l|) (4.7) 

must be satisfied in order for acceleration to occur. The loss-cone angle is defined by aic = 
sin"' [(Sup/5down)''^^], where Sup and Sdown denote the magnitude of the magnetic field in 
the up- and downstream region, respectively. 



Px4 




Fig. 4.4 Illustration of shock-drift acceleration in the — plane for a shock with ft = 0.8 and a/^, = 45°. 
All particles initially located in the light grey area are transformed into the dark one. 



Shock-drift acceleration represents a transformation in the jSx — jS|| plane as illustrated 
in Figure l4~4l As seen in Figure l4~4l the accelerated distribution function is simply a shifted 
loss-cone distribution. The flux of the accelerated electrons parallel to the upstream mag- 
netic field is defined by 

=2;r^oc4 j^'d/3|j/3|| jJ^'yrf/3x/3x-/o..(/3||,/3x), (4.8) 
and the differential flux is given by 

dFaccW 

y'flcc, II ~dE ' (4.9) 

Figure l43] shows the differential fluxes obtained from the shock-drift acceleration model 
for various angles d, assuming an initial distribution in the form of a Maxwellian with a 
temperature 7" = 10 MK. For comparison, the dotted line in Figure |43] represents a pure 
Maxwellian distribution with the same temperature. One can see that shock-drift acceler- 
ation indeed provides a significantly enhanced component of energetic electrons, particu- 
larly for large values of the angle d. These analytically-obtained results are qualitatively 

confirmed by num erical simulations performed for plasma conditions usually found at the 

Earth's bow shock jKrauss-Varban & Wull989i : lKrauss-Varban et al .ll989l : lKrauss-Varban & Burgesd 
[1991). 

In order to compare these results with observations, some basic paramet ers must be 
provided. Here, SOL2003-10-28T11:10 (X17.2) is considered as an example ( iMann et"al] 
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Fig. 4.5 Differential flux of electrons accelerated by shock-drift acceleration for various shock angles. The 
dotted line represents the flux of a pure Maxwellian distribution. 



l2009l) . In this event, the termination shock appears at 300 MHz for the harmonic emis- 
sion. Matching harmonic emission to the observed 300 MHz frequency requires an electron 
number density of 2.8 x 10** cm^-*. According to a two-fold density model of the corona 
( lNewkirlJll96lh . such a region is located ^160 Mm ab ove the photosphere, at which a mag- 
netic field of ~5 G is expected & McLeanll 19781) . These parameters lead to an Alfven 
speed w 560 km s^'. For the temperature, a value of 40 MK is adopted, resulting in a 
thermal electron speed of 2.5 x 10^ km s^' . 

Assuming a density jump of a factor of two across the termination shock, the Rankine- 
Hugoniot relations give a magnetic field S(jo„n/Sup = 2 and an Alfven-Mach number Ma 
= 3.5, corresponding to a shock speed of ~2000 km s^'. Due to the curvature of the 
magnetic field lines convected through the termination shock, there is a distribution of 9 
across the shock. Recognizing from Figure l43] that the flux is very sensitive to 8, we model 
this curvature by a circle; the resulting flux of accelerated electrons can be integrated over 
6, giving an accelerated electron production rate = 2 x lO-'^ s^' and a related power 
Pf = 1 X 10^^ erg s~^ above 20 keV. This agrees well with the results inferred from RHESSI 
observations dWarmuth et al.ll2007l,l2009l) . 

It must be stressed that a non-negligible part of the energy of the outflow jet is trans- 
ferred into accelerated electrons. Thus, such a test-particle approach has its limitations: the 
termination shock acceleration process must be treated in a fully kinetic manner, in which 
the accelerated particles themselves are included in the conservation laws across the shock. 



4.4 Particle acceleration in a single 3-D RCS with complicated magnetic topology 
4.4.] Basic equations 



Following Section [33l we consider a highly-stressed current sheet oriented in the x-direction, 
with a lateral field 

B,{x)=BQtanh(^-^^ (4.10) 
and a t ransverse magnetic field component that is dependent on z according to dzharkova & Gordovskvvl 



B,{z)=Bo(-] , (4.11) 
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where a > (we consider an illustrative value a = 1). A weak, constant longitudinal com- 
ponent of the magnetic field By is also assumed to be present: 

By = pBo. (4.12) 

Particle trajectories are calculated from the relativistic equations of motion 

dr p 

dt ihqJ 



(4.13) 



"^P "'E+-— xB), (4.14) 



dt \ c moY 

where t is time, r and p are the particle position and momentum vectors, E and B are the 
electric and magnetic field vectors, q and mo are the charge and rest mass of the particle, c is 
the speed of light, and y = \/ (p/moc)^ + 1 is the Lorentz factor. The integration time steps 
are chosen to be much smaller than the corresponding gyro-period, i.e., <^f < 0. 1 (in/q) B^^ . 
For protons, this typically requires a time step < 10^^ s, while for electrons the required 
time step is substantially shorter, <4 x 10^" s. 

It should be noted that in the above description, the electric and magnetic field vectors E 
and B are prescribed quantities; the self-consistent electric and magnetic fields produced by 
the accelerated particles thems elves are neglected. Such a test particle approach is valid for 
up to only ~10^ test particles ( iBirdsall & Langdoiilll985h . and we shall return to this point 
in Section 14.71 It should also again be noted (see Section 13.31 1 that in this simulation the 
reconnection electric field is assumed to extend over the entire region covering the whole 
diffusion region. Thus, the caveats noted in this earlier Section also apply here. 

4.4.2 Particle motion inside a 3-D RCS 

In the test-particle description used here, the electron trajectory and the direction of ejection 
are determined uniquely by its initial coordinate, velocity and the magnetic field geometry. 
Particles inflowing into a 3-D RCS are accelerated mostly near the lateral field reversal at 
the jr-axis. In the presence of a non-zero By, electrons and protons are ejected in opposite 
directions from the midplane (x = 0); whether the particle is ejected to ;c > or .v < 
depends on the sign of the particle charge, and the signs and magnitudes of B^ and By. If 
By is strong enough, then all protons (regardless of the side of the RCS they entered from) 
are ejected to the x> semispace, while all electrons are ejected to the opposite semispace, 
as illustrated in Figure l4~6l ( Zharkova & Agapitov 2009). For particles of the same charge, 
there are two fundamental types of trajectories inside an RCS: particles that enter from and 
then are ejected into the same semispace are hereafter referred to as "bounced" particles, 
whereas particles that are ejected into the opposite semispace are referred to as "transit" 
particles. 

Electron/proton asymmetry in the particle traje ctories is dependent on the magnitude 
of the guiding field jzharkova & Gordovskvvll2004l) . The asymmetry is full (equal to unity) 
if By is strong enough (5y/S, > 1.5 x 10^^) that electrons are ejected into one semiplane 
and protons to the opposite one. If the guiding fie ld is weak {B^,/B^< 10~^), electrons an d 



protons can be ejected into the same semiplane jZharkova & Gordovskvv 20041 . l2005bl) . 



For a guiding field of intermediate magnitude, electrons and protons are ejected in a par- 
tially neutralized manner, with protons dominating in one semiplane and electrons in the 
other Izharkova & Gordovskv vl f2005bl) . The result of this separation of particle species is 
a polarization electric field appearing between the protons at the edge and electrons at the 
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Fig. 4.6 Proton (blue lines) and electron trajectories (purple lines) in a 3-D RCS w i th a la teral magnetic field 
B.,- = 10 G and an drift electric field Ey = 100 V m"' . From lZharkova & AgaDitwl j2009h . 



midplane of the RCS. Suc h a field was anticipated bv iLitvinenko & Somov (Il993ah and 
lOreshina & Somovl ( l2000bh . and has been investigated in some detail bv lzharkova & Agapitovl 
( l2009l) . 

The model produces power-law energy distributions with spectral indices 3 = 1.5 for 
protons and 5 = 2 for electrons. Such a result is confirmed by both analytic i Litvinenkq 
ll996cJ) and numerical approaches dzharkova & Gordovskv\l2005bl :IWood & Neukirch' 20051). 
For com parison, the particle energy spectra calculated for the model RCS of Oreshina & Somovl 
( l2000bl) and presented in Figure l4~7l show features beyond such simple power-law forms: 
there is a rapid rise in flux to a maximum at £ = above this, the spectrum is approx- 
imately power-law. £„, 10 — 12 keV for protons (ci 2 — 5 keV for electrons); these val- 
ues increase slightly with the value of the parameter jS = B^./B,. The energy E„, can be 
considered as the lower energy cutoff to the electron spectrum (see Section 12.11 1. and the 
values obtained here agree quite well with the values deduced from RHESSI observations 
jHolmanetal .11201 ih . 

The spectral index 5 for the power-law part of the spectrum (E > £„,) is found to vary 
slightly with the magnitude of the guiding field By. For a weak guidi ng field ([5 = 10~"* - Fig- 
ures l4.7l [a] and [d]), the particle trajectories are nearly symmetrical dzharkova & Gordovskvvl 
i0O5b); the particles are ejected equally into both semispaces as approximately neutral 
beams. The particle energy spectra in both semispaces (x < and ;c > 0) are very simi- 
lar, with a sharp increase from zero to £„,, followed by a power-law with spectral indices 
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Fig. 4.7 Energy spectra of electrons (dashed lines) and protons (solid lines) ejected from the RCS. Spectra (a) 
- (c) are calculated for /5(= By/Bo) equal to lO"'*, 10^^ and 1, respectively, and correspond to the semispace 
X > 0. Spectra (d) - (f) are calculated for the same values of /3, but for the semispace x < 0. The estimated 
proton spectrum is presented as a dotted line in plot (c) and the estimated electron spectrum is shown as a 
dash-dotted Une in plot (f). From .Zharkova & Gordovskvv, t2005h) . 



5 1.8 for protons and 5 ~ 2.2 for electrons. The power-law shapes extend up to about 
100 keV (electrons) or 1 MeV (protons). 

For a moderate guide field (jS = 10^^ - Figures l477l rb1 and [e]), the protons and electrons 
are also ejected into both semispaces; however, the symmetry of their trajectories is partially 
destroyed. Hence, most protons are ejected into the semispace x > with 5 1.7, while a 
few low-energy protons are ejected into the semispace x < with very soft, thermal-like, 
spectra (5 ~ 4.8). The opposite picture holds for electrons: the bulk are ejected into the 
semispace jc < with an energy spectrum 5 ~ 2.0, while a smaller number of the low- 
energy electrons are ejected into the semispace x> with a thermal-like spectrum having a 
maximum £,„ w 1 keV. 

Finally, if the guiding field is very strong (/3 ~ 1 - Figures l477l [c] and [f]), the particle 
trajectories are completely asymmetric. Particles are ejected completely separately from the 
RCS midplane into opposite semispaces: all protons are ejected into the semispace jc > 
while all electrons are ejected into the semispace x < 0. Their energy distributions show a 
sharp increase of particle number up to around a few keV followed by a power-law energy 
spectrum of index 5 :^ 1 .5 for protons and 5 :^ 1 . 8 for electrons. (These spectral indices vary 
with the magnitude of a guiding field (Fig. 14.8b reflecting different reconnection scenarios 
and with different phases of magnetic reconnection, as discussed in Section l43] ) 

As shown in Figure 14.91 during the reconnection process, the shape of the reconnect- 
ing current sheet changes from a non-compressed one with a symmetric X-type null point 
to a strongly-compressed shape with a well-defined elongation along the z-direction and a 
compression in the jc-direction. The extent of this deformation depends on the value of the 
parameter a jPriest & Forbesll200C Sornovll2000b and possibly corresponds to different re- 
connection rates ( Arber & HavnesI 20061) . At the reconnection onset (a = 0) the diffusion 
region is located very near the null point. As the reconnection progresses (a increases from 
to 1), the diffusion region becomes much wide r in both x and z-directions. A s a increases 
further, the reconnection rate begins to decrease dZharkova & Agapitovll200^ . 
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Fig. 4.8 The energy spectra resulting from the protons entering the cuirent sheet from opposite sides, fo r 
different ratios j3 = By/Bo and different magnitudes of the parameter a. From lZharkova & Agapito^ )2009l) . 



The spectral index of the accelerated particle energy spectra is associated with the 
parameter a and depends on the magnitude of the guiding field (Fig. 14.8b . For an RCS 
with a density that increases exponentially wit h distance z and for a variat ion of density 
n{z) °^ {z/a)^ , the following results are obtained dZharkova & Agapitovll2009h : for electrons 
in a weak (strong) guiding field By, 



C^e.weak " 



1 + - 



2 + X ^ 
la 



1 1+A 

e>.,strong - 2 + 2a ■ 



while for protons in a guiding field of any strength, 



1 + - 



1+A 
la 



(4.15) 



(4.16) 



In principle, therefore, by simultaneously measuring the spectral indices of electrons and 
protons in the same flare, one can deduce the magrietic field geometry under which th ese 
particles were accelerated teharkova & Gordovskvvl l2005bl : Izharkova & Agapitovll200^ . 

We again stress that all of these results have been arrived at using a test-particle ap- 
proach, in which the feedback of the electromagnetic fields associated with the accelerated 
particles is not taken into account. Such a limitation is addressed in the PIC approach dis- 
cussed in Sectionl4~8] 



4.5 Particle acceleration in 3-D MHD models with fan and spine reconnection 



A further development in modeling particle acceleration in a reco nnecting current shee t 
involves more advanced MHD models of fan and spine reconnection jPriest & Titov|[l99^ . 

The axis of symmetry of the magnetic field, called the spine and designated here as the 
X axis, is a critical field line connecting to the null point. In the plane x = 0, the magnetic 
field lines are straight lines through the null point describing a fan; hence, the (vjz) -plane is 
termed the /an plane. The spine and the fan are the 3-D analogs of the separatrix planes in 
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Fig. 4.9 Views of the current sheet for differen t values of a = 0.5; 0.7; 1.0; 1 .5; 2.0 and 5.0 (plotted by rows 
from the top left to the right bottom plot). From lZharkova & Agapitovl i2009h . 



the classic 2-D X-point geometry. The field configuration near the potential magnetic nulls 
is given by Eauations l4. 101 - 14. 121 

Two regimes of re connection at a 3-D magnetic null are possible: spine re connection 
and fan reconnection l lLau&Finrj[l990l ; IPriest & Tito^[l996) . While the magnetic field 
configuration is the same in the two regimes, they are characterized by different plasma flow 
patterns and different electric fields. Spine reconnection has a current concentration along 
the critical spine field line, while fan reconnection has a current sheet in the fan plane. 

For a potential null, the spine reconnection regime is characterized by plasma flows lying 
in planes containing the spine, while fan reconnection has non-planar flows that carry the 
magnetic field lines in a "swirling"-like motion symmetric about the fan plane. It should be 
noted that the electric field and flow pattern are non-axisymmetric in all cases. The model of 
|Priest&Tito"vl ( ll996h provides a useful analytical model for initial studies; it allows some of 
the key features of topology and geometry to be identified, and it is the closest 3-D analog of 
the widely-studied constant out-of-plane electric field used in 2-D studies. It describes only 
the outer ideal reconnection region and does not include acceleration due to the parallel 
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Fig. 4.10 Traject ories and energy gains in f an reconnection, for particles injected on tlie boundary at different 
latitudes /3. From lDalla & Brownind <2008l) . 



electric field in the inner dissipation or resistive region; indeed, it breaks down very close 
to the spine/fan, as the electric fields formally diverge. However, since few particles will 
enter the small-volume dissipation regions, the results are likely to be appropriate for the 
bulk of the energy spectrum. This model has been use d as a basis for particle acceler ation 
at 3-D nulls in th e regimes of both spine r econnection jPalla & Brownin3l2005l l2006h and 
fan reconnection jPalla & Brow ning 2008). Particle trajectories are obtained numerically by 
solving the relativistic equations of motion (Equation l4.14b . Again, a test-particle approach 
is used (for limitations, see Section|4j7}. 

The strength of the driving electric field is quantified by the parameter 

A = ^, (4.17) 

where vx is the perpendicular speed associated with the gyro-motion, which may be equated 
to the thermal speed of the plasma. This can also be expr essed as the dimensionless m agnetic 
moment when speeds are normalized to the drift speed dVekstein & Brownindll997h : 

inv^, /2Bo 

fi= z^n/.R ' (4.18) 
/w| (L)/2Bq 

where Bq is the amplitude of the magnetic field. If /i <S 1 , the electric drift speed at global 
length scales is strong compared with the thermal gyro-motion and we are in the strong drift 
regime corresponding to fast reconnection; in practice, this regime is reached when p L w 1. 
In this case, particles may undergo significant acceleration jPalla & Browrun3l2005h even 
in the outer "ideal" reconnection region, similar to those discussed in Section l4!4l 

Typical individual trajectories in the fan reconnec tion are shown in Figure [ 4. 101 while 
similar plots for spine reconnection can be found in iDalla & Brownin j ( |2005|) . Note the 
strong dependence of the energy gain on the injection location, similar to those found for a 
single RCS in Section l4~4l particles injected closer to the fan plane (at lower latitudes jS) are 
accelerated the most as they approach the strong-electric-field region near the fan plane. 

Energy spectra for the spine and fan reconnection (Figures [4.1 ll and [4.12l respectively) 
show that a significant fraction of protons can be accelerated to energies of up to 10 MeV, 
comparable to the values obtained in the simple 3-D topology presented in Section 14.41 
As appears to be natural during any reconnection process, a generally power-law spectrum 
is obtained for the accelerated particles, but with a significant "bump" at higher energies. 
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Fig. 4.11 Energy spectra of particles in spin e reconnection at successive times, showing evolution towards 
steady state. From lDalla & Brownind )200d) . 
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Fig. 4.12 Energy spectra of particles in fan reconnection at successive times, showing evolution towards 
steady state. From .Dalla & Browning 1 2008.) . 

particularly for the spine reconnection, for which the bump arises from particles approaching 
close to the null itself. 

This "bump-on-tail" distribution is similar to those reported for a single 3-D current 
sheet in Figure l4~8] in Section l4~4l (see also lZharkova & Gordovskvvl2005cl : lzharkova & Agapitovl 

and is formed by the two populations of protons: "transit" and "bounced" ones, re- 
spectively. Note that the slopes of the spectra are similar for the spine and fan reconnection 
in the intermediate energy range, although fan reconnection is significantly less effective at 
generating high-energy particles. 

Detailed information on both the spatial distributions of the accelerated of particles, 
and where these particles originate, are shown in Figures [4.131 and [4. 141 for spine and fan 
reconnection, respectively. These show the angular location of the particles at the initial (top 
plot) and the final (bottom plot) times. Each location is identified by its latitude jS, i.e., the 
elevation angle from the x-y plane, and longitude (j). Inflow regions are the top-right and 
bottom-left quadrants in the (<^, jS) representation, but for clarity, only particles originating 
in the latter are shown (the other quadrant can be added by symmetry). Each point represents 
one particle and is color-coded according to the final particle energy. 

In the case of spine reconnection, particles injected in the plane in which the electric 
field is weakest (longitude (j) = 0) actually have the gr eatest energy gain, although their 
acceleration times are longest ( iDalla & Browninal2005h . The high-energy particles which 
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Fig. 4.13 Positions of particles at initial and fi nal times for spine reconnection, color-coded according to their 
final energies (from lDalla & BrowninjlaOOd) . 
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Fig. 4.14 Positions of particles at initia l and fi nal times for fan reconnection, color-coded according to their 
final energies. From lDalla & Brownind <20O8l) . 



escape mostly emerge in jets along the spine. However, most of the highest-energy particles 
actually remain trapped in the magnetic field (Dalla & Browning 2006). The occurrence of 
these high-energy tra pped particles may be an explanation for coronal hard X-ray sources 
dKrucker et al]|2008ah . 

In fan reconnection, particles starting near the fan plane are those that gain the largest 
energy, as would be expected since the fan plane is the location where the electric field 
is largest. As one moves towards increasing latitudes, particles will reach the fan plane 
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with greater difficulty and so gain less energy. Figure 14.141 shows that, at the final time, 
particles are confined to a range of longitudes [-130°, -50°]. This is a result of the tendency 
of the electric field drift in fan reconnection to push particles towards the x = plane after 
they have passed the region near the null. Figure l4~T4l also shows that the initial azimuthal 
locations near (j)=90° are the most favorable for particle acceleration since, near this plane, 
the transverse component of the electric field drift, pushing particles towards the fan plane, 
is largest. The fan reconnection regime is less efficient in seeding test particles to the region 
of strong electric field, compared with spine reconnection, and acceleration times are longer. 

These results, obtained for a single-current-sheet recon nection model with a m ore com- 
plicated magnetic and electric field topology simulated by 'Priest & Titovl \l99&), advance 
further the conclusions derived from the basic 3-D reconnection models discussed in Sec- 
tion l4.4l This approach confirms that acceleration by a super-Dreicer electric field generated 
during magnetic reconnection is a viable mechanism, which naturally explains the energy 
input from magnetic field dissipation and its output in the form of power-law energy spectra 
of accelerated particles, with spectral indices related to the magnetic field geometry occur- 
ring during the reconnection. Apparently, this approach offers a good diagnostic tool for the 
magnetic field topology from the observed signatures of accelerated particles. We yet again 
note that the limitations of the test-particle approach (see Section l477] ) can restrict this diag- 
nostic capability, and that a full kinetic approach to particle acceleration in such geometries 
is required. 

4.6 Particle acceleration in multiple current sheets formed by MHD turbulence 

Although particle acceleration in a single current sheet can, in principle, account for some 
flares or even flare ribbons, it is highly likely that flares are due to a number of unresolved 
magnetic tubes interacting with each other, so that more complex coronal geometries with 
multiple acceleration sites must be considered. One obvious benefit of such a scenario is that 
a particle can be accelerated at more than one location, perhaps enhancing the energization 
process. 

4.6.1 Acceleration in Cellular Automata Models 

The idea of multiple dis sipation sites in flares is longstanding, but the first serious effort 
at modeling this was bvl Lu&Hamiito]l (fl99lh. who took the ideas associated with self- 
organized criticality (SOC) and applied them to the solar corona. They showed that if the 
right conditions were met, the triggering of dissipation at a single point could lead to a 
"spreading" of dissipation across a large coronal volume and that the distribution of event 
size as a funct ion of energy follow ed a power law {■yK{S') °c S'^^'^) similar to that observed 
in flares (e.g.. lCrosbv et alj|l993[) . However, a fundamental question remains: can existing 
MHD models verify that the main rationale behind the SOC theory is valid? 

Self-organized criticality models identify diss ipation occurring at many spatially-separated 
regions as illustrated in Figure I4TT5] taken from lviahos et all ( l2004h . The upper left panel 
shows a coronal field geometry, reconstructed from photospheric magnetograms, in a vol- 
ume with characteristic dimensions of 10^ cm. The current density in a localized region is 
shown in the upper right panel, and a snapshot of the locations where currents exceed the 
threshold for dissipation in the cellular automaton model (referred to as Unstable Current 
Sheets) is shown in the bottom right panel. Energy release in this model occurs at a number 
of locations of unstable current sheets (bottom left panel). 




Fig. 4.16 Particle trajectories in tlie CA model. From lVlahos et alj i2004) . 
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To assess particle acceleration in such a coronal geometry, Iviahos et all ( l2004h have 
developed a model in which the particles move from one unstable current sheet to another 
unstable current sheet, gaining (or losing) energy at each and undergoing ballistic motion in 
a background magnetic field in between. This is illustrated in Figure l4T6l where a particle 
first interacts with a site escapes along the magnetic field, encounters a site /+ 1, escapes 
again, and so on. The energy gain thus occurs in many small increments, rather than at the 
single location of the current-sheet approach. By using the self-organized criticality hypoth- 
esis, one can identify the structure of the spatial distribution of the unstabl e current sheets 
inside a complex active region which exhibits a specific fractal structure jMcIntosh et alJ 
l2002h . Knowledge of the spatial structure of the unstable current sheets can help us to recon- 
struct the probability distribution of its spatial scales, and th us, the probability of a particle 
colliding with a given element of it ( ilsliker & Vlahosll2003h . We do not yet understand the 
statistical properties of the "kernel" - the energy gain or loss for a particle colliding with 
an unstable current sheet - but this is relatively easy to study using three-dimensional MHD 
codes. 

We can therefore start from a global reconstruction of the distribution of unstable cur- 
rent sheets and recover the kinetic properties of the particles, under the assumptions of the 
self-organized criticality theory. : The acceleration is thus determined by three probability 
densities: one, defining the distance between a pair of accelerators; the second, the electric 
field strength at each accelerator; and a third, either the time spent in, or the length of, the 
unstable current sheet. The first two are assumed to be power laws, and the third is taken 
to be a Gaussian. Note that this model only addresses acceleration by direct electric fields, 
since that can be simply parameterized at the dissipation sites. 

In the example shown in Figure |4~T5] the particle moves in a box of size lO''* cm and 
the distance between accelerators ranges from 10"* to 10'" cm. The magnetic fields may be 
such that the particles become partially trapped and the trajectory becomes rather complex. 
The upper limit of the free travel distance is set to a value larger than the coronal volume 
in order to allow particles to leave the coronal region without undergoing any interaction 
with electric fields, besides the initial one. Note that for the maximum time At = I second 
for which the system is monitored, these escaping particles will move at most a distance 
cAt = 3 X 10'" cm, which is of the order of the length of the coronal volume considered. 

The electric field at each unstable current sheet varies between sub-Dreicer and 10^ 
times Dreicer, and the acceleration time has a mean value of 2 x lO^"* s ( Vlahos et al.. 2004h . 
A large number of test particles are tracked through this coronal geometry using the rela- 
tivistic equation of motion described by Equations 14. 141 with acceleration taking place very 
rapidly (<C 1 s). 

Figure 14.171 shows the distribution function of electrons at four different times. The 
initial distribution is a Maxwellian with a temperature of 100 eV (dash-dot line). One can 
see that (a) prompt acceleration to energies well in excess of 100 keV takes place, (b) just 
under half of the particles end up in the high energy tail and (c) the tail is a power law with 
slope 5 ~ 4. When the test particles are protons, much longer acceleration times must be 
considered to allow for similar energies to be reached within the prescribed acceleration 
time. 

Another approach based on the SOC models was proposed bv lAnastasiadis et"al] ( |2004|) 
to predict the kinetic energy distributions of electrons. The acceleration of electrons is still 
due to the presence of the randomly placed localized electric fields produced by the energy 
release process and simulated by the CA model. It is assumed that electrons interact succes- 
sively with reconnecting current sheets with the electric field linked to the energy release 
time series. At each encounter, they randomly gain or lose energy. Electrons can be acceler- 
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Fig. 4.17 Accelerated electron spectra in the CA model jVlahos et al.l2004t) . 

ated from a Maxwellian initial distribution up to a few MeV, and can be represented either 
by a power-law or an exponential distribution. 

As the number of interactions increases, the electron energy distribution begins to di- 
verge from a well define d power law, and an exponential tail develops (e.g., iTrottet et al.l 
1 19981 : IVilmeretal.ll2003h . Thick-target radiation from these energetic electrons was also 
computed, leading to hard X-ray spectra at low energies (10-100 keV) which are much 
flatter than usually observed. At higher energies (100-1000 keV), the X-ray spectra were 
consistent with the observations at these energies. 

4.6.2 Acceleration in three-dimensional MHD turbulence 

Numerical experiments for test-particle acceleration in the turbulent magnetic and elec- 
tric fields obtained from pseudo-spectral direct numerical solutions of the compressible 
three-dimensional M HD equations, with the Ha l l term included, have been carried out by 
iDmitruk et all ( |2003|) and lDmitruk & MatthaeusI 1200^ . In both experiments electrons and 
protons were accelerated within very short times before particles left the simulation box 
(with length size of the order of a few turbulent correlation lengths). The momentum distri- 
bution functions very quickly develop long high-energy tails. 

For electrons, there is no substantial difference found between the results with and with- 
out the Hall term in the MHD solution, while for protons significant differences are seen 
in the low-momentum part of the distributions. When a background uniform magnetic field 
(the guiding field) is added, the particle acceleration is anisotropic. Electrons develop large 
parallel momentum, while protons develop large perpendicular momentum similar to those 
reported for a single current sheet in Section l4~4l 
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4.6.3 Acceleration in stochastic current sheets 

iTurkmani etal] ( |2005| , |200^ have i nvestigated the issue of partic le acceleration in a 3-D 
coronal MHD model developed by iGalsgaard & Nordlund il996l) . The initial setting for 
the model is a straight stratified coronal loop with a magnetic field which extends between 
the two photospheric footpoints. The response of the loop to random photospheric motions 
is then followed. A turbulent cascade develops, leading to a formation of multiple current 
sheets in the corona. The resulting electric fields comprise both the inductive electric field 
and the resistive one, the lat ter appearing from the hybrid resistivity used in the model of 
IGalsgaard & Roussevl feOQ21. However, it is found that the impact of the (mainly perpendic- 
ular) inductive field on acceleration is weak. The distribution of current sheets changes in 
time with the turbulence development, giving rise to the possibility of very bursty particle 
acceleration events. 

In order to study particle acceleration in this enviroimient, the magnetic field topology 
deduced from the MHD simulation is assumed to be quasi-stationary, similar to that in Sec- 
tion [33] Then, test particles are followed through their motion in the electromagnetic fields 
associated with the multiple current sheets, using the relativistic equations of motion (Equa- 
tion [TTlll. The model also allows particles to leave the loop once they hit the boundaries, 
thus limiting their energy gain. The simulation results show different types of particle dy- 
namics, with particles both gaining and losing energy; however, not surprisingly, most leave 
the loop with net energy gain. About one quarter of the particles do not interact with the 
current sheet at all and thus gain no energy. 

The results show that the resistive electric field, being a parallel field, is the princi- 
pal means of accelerating particles. The acceleration is very fast and efficient, accelerating 
particles to relativistic energies in a very short time. The resulting distribution function of 
accelerated p articles has a clear do uble-power-law structure, with a thermal component at 
low energies ( iTurkmani etani2006l) . 

The MHD model currently in use produces particles with higher ener gies at some- 
what shorter times compared to the Cellular Automata model discussed above dviahos et alj 
l2(K)4h . This is likely because of different choices in the basic parameters, especially the 
magnetic field strength and density. The same MHD model accelerates both electrons and 
protons without any need to adjust the electric fields, contrary to the CA-based model. The 
double power law does not arise in the CA model because of the possibility of particle loss 
in the MHD model. Both approaches show that fragmented current sheets in the corona are 
a reasonable scenario for particle acceleration. 

iDecamp & Malaral ||200^ performed a test-particle numerical experiment to simulate 
particle acceleration in the low-frequency turbulent environment generated by footpoint mo- 
tions in a coronal loop. Only the effect of the resistive electric field E is retained, which is 
mainly parallel to the axial magnetic field. In the spatial frequency spectrum, the contri- 
bution of small scales is dominant. The spatial structure of E is obtained by a synthetic 
turbulence method (p-model), which allows intermittency with a power-law distribution to 
be reproduced. Spatial intermittency plays a key role in acceleration, enhancing both an 
extension of a power-law range, and an increase in the maximum energy. 

By solving the relativistic equations of motion, the time evolution of the particle distri- 
bution can be determined. Electrons are accelerated to energies of the order of 50 keV in 
less than 0.3 s, and the final energy distribution also exhibits a power law. A correlation is 
found between heating events and particle acceleration that is qualitatively similar to that in 
the stochastic models discussed in Section [34l 
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4.6.4 Limitations ofMHD turbulence models 

With all their attractiveness in accounting for observed electron distributions in solar flares, 
models invoking MHD turbulence in an ensemble of current sheets nevertheless have a few 
essential problems. 

First, MHD models involving multiple current sheets inevitably inherit the same size- 
and time-step problems as the PIC simulations discussed in Section 14.81 allowing us to 
simulate MHD turbulence and electron acceleration in only a very small region (about 1 m-'). 
It is then assumed, without justification, that the results from this small volume can be 
extended to the coronal scale. 

Second, the basic model of a current sheet (discussed in Section 13. 3t is still used for 
the energy gains by accelerated particles in these reconnecting current sheets, while more 
realistic 3-D models of the reconnecting current sheets (as discussed in Sections r4.4l and l431 l 
show the energy gains to be significantly different. 

Third, even if, as suggested above, multiple current sheets are formed in the whole flare 
volume, the diffusion regions in each have to be connected somewhere to reconnecting mag- 
netic field lines. In a single current sheet model (see Sections r3.3l and l4~4l l. and as also shown 
in the origin al magnetic field restoration from the magnetograms taken at the photosphere 
dviahos et al. 2004). these lines are embedded in the photosphere. The trajectory simula- 
tions presented in Section l44l show that the particle drifts follow the magnetic field lines of 
the main magnetic component and move towards the photosphere, not parallel to the solar 
surface. Therefore, in the proposed multi-current sheet models, the issue arises as to how the 
test particle can mo ve across the magnetic fi eld, in order to get to the another current sheet 
(see, in this context, lEmslie & Henouxlll995h . 

Last, but not least, a major problem arises from the fact that these simulations inherit the 
most common limitations related to the test-particle approach, as discussed in Section l4~7l 



4.7 Limitations of the test particle approach 



As discussed above in Sections 13.21 (basic reconnection in a single 2-D current sheet), 13.41 
(acceleration by shocks and/or plasma turbulence), 14.41 and 14.51 (reconnection in 3-D ge- 
ometries), and l4.6l (MHD models incorporating multiple reconnection sites), a considerable 
amount of simulation of the particle acceleration process involves a use of "test particles" 
moving in a prescribed electromagnetic field configuration. Such a test-particle approach is 
certainly useful at some level - for example, it allows us t o investigate the individual tra- 
jectories of electrons and protons inside 3-D curre nt sheets Jzharkova & Gordovskvvll2004l : 
iDalla & Brownindl2005l ; lwood & NeukiTd]|2005h. and to obta in the aggregat ed energy dis 



tribut i ons for a large volume of particles dZharkova & Gordovs kvv 2005a, b; D alla & Brownind 



l2005l : lBrowning & DalliJl2007h . which can be compared with those inferred from observa- 
tions. The approach can be used iteratively to extend the results of exact PIC simulations to 
a much larger region (see Section l4~8l l. 

However, considerable caution must be exercised when applying a test-particle approach 
to situations in which the number of accelerated particles is sufficiently large (see Sec- 
tion [ZT3]l. In such cases, the electric and magnetic fields associated with the accelerated 
particles themselves rapidly become comparable to the fields pre-supposed at the outset. 
Recall (Section lZTT l that a major flare can accelerate some 10^^ electrons s^' and that these 
electrons propagate within a magnetic flux tube of order lO' cm in radius. Attempting to 
model this as a single large-scale ejection leads to the well-known (e.g., lHolman|[l98l ; 



Particle Acceleration 



47 



lEmslie&Henouxlll995h paradoxes that (1) the magnetic fields induced by the resulting cur- 
rent vastly exceed the pre-supposed field, and (2) that the electric fields resulting from both 
charge separation and induction effects associated with the accelerated particle stream (see 
Sections [4.4l and l4~8l l exceed plausible values by even more absurd margins. Specifically: 

- Ampere's law, applied to a steady flow of lO-'^ electrons s^' (current I ^ lO'^ Amps), 
in a cylinder of cross-sectional radius r w 10^ m gives a corresponding magnetic field 
strength B w fioI/2Kr w 3 x 10^ T w 3 x 10** Gauss. Using a propagation volume V w 
10^' m"*, the associated energy content (B^ /2y.o)V w lO"** J w 10^^ ergs, at least ten 
orders of magnitude greater than the energy in the electron beam itself; 

- The typical rise time for a hard X-ray burst in a solar flare isTwl — 10s. To initiate a 
current / over a timescale T in a volume of inductance L w p.ol (where £ is the charac- 
teristic scale) requires a voltage V w /Xo£7/t. Alternatively, the largest current that can 
appear when a voltage V is applied is 7 w Vt/ figi. With the characteristic voltage asso- 
ciated with the acceleration of hard-X-ray-producing electrons being ~30 kV, and with 
solar loop sizes £ w 10^ m, this gives a maximum current / w 10^ Amps, some fourteen 
orders of magnitude less than the total current involved. 

It is clear, therefore, that the self-consistent electromagnetic fields associated with electron 
acceleration rapidly become comparable to the fields pre-supposed to be present in the ac- 
celeration region; in such an environment, a test-particle approach breaks down severely. 
Given this, it is a clear imperative to take into account the self-consistent evolution of the 
electric and magnetic fields, including the fields associated with the accelerated particles 
themselves. Failure to do so can result in unrealistic numbers of accelerated particles pre- 
dicted by acceleration models and even paradoxical results, such as the energy in accelerated 
particles exceeding the free energy in the pre-flare magnetic field configuration. 

4.7.1 The polarization electric field 

As discussed in Section 14.41 (see Figure I4.6I I. reconnection in a 3-D current sheet geome- 
try leads to a separation of protons and electrons across the midplane. As shown in Figure 
14.181 such a charge separation results in a substantial (up to 100 V cm^') Hall-type polar- 
ization electric field across the midplane. As pointed out in Section [T2l inclusion of such 
Hall-type electric fields can significan tly increase the effectiveness of the whole magnetic 
reconnection process tBirn et al..,200ll ; lHuba & Rudakovll2004l ; lDrake et al.ll2005h ; see also 
Sectioning] 

Electrons, because of their much smaller gyroradii, typically do not show much change 
in behavior in the presence of the polarization field. On the other hand, protons, with their 
larger gyroradii, gain slightly more energy when they move in the positive jc-direction than 
in the negative jc-direction. The net (nonlinear) result is that protons gain energy as a result 
of the polarization field, thus tending to make the proton energy spectrum flatter. A full 
investigation of the effect of the polarization electric field requires use of a PIC approach, as 
described in Section l4~8] 

4.7.2 Turbulent electric fields 



Various turbulent electric fields typically form inside a reconne cting current sheet - these 
are caused, for example, by the two-stream inst ability ( Bretll2009l) of the two electron beams 
that enter the current sheet from opposite sides jZharkova & Agapitovll200^ . Such electric 
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Fig. 4.18 Polarization electric field induced by the separation of protons an d electrons accelerated in a 3-D 
current slieet with Bo= 10 G, at various distances z, from the null point ffrom lZharkova & Agapitovll2009l) . 



fields, and their associated currents, typically oscillate near the plasma frequency and have 
an exponential growth rate F and an energy saturatio n level PE given by the analytical 
relations derived for the cold-beam plasma instability jShapiro & Shevchenkol Il988l ; iBrell 
12009) : 



/T73 2-4/3 4^/3 



CO, 



'pe-i 



PE = {A/2fl^Eb,, 



(4.19) 



where A is the ratio of the beam to plasma densities, and Zibeam is the beam kinetic energy. 
The effect of the electric fields associated with such wa ves on the accelerated particle distri- 
butions was estimated bv lZharkova & Agapitovl l l2009l) : this work confirmed the estimates 
of li^boogl) . Further progress requires full kinetic simulations. 



4.8 Particle-in-cell simulation of acceleration in a 3-D RCS 

As has been pointed out repeatedly above, a full description of the particle acceleration 
process in solar flares requires that the electromagnetic fields and currents associated with 
the accelerated particles be self-consistently included. In this section, we describe early 
efforts to accommodate this crucial element of the modeling. 



4.8.1 Description of the PIC approach 

In order to self-consistently include the electric and magnetic fiel ds induced by the ac- 
celera ted particles, one can use the 2D3V (PIC) simulation code jVerboncoeur & Gladdl 
Il995h . This PIC method solves the equations of motion for the plasma particles (see Equa- 
tions |4T4]|, using values for the electric and magnetic fields which include both the pre- 
specified background fields E and B and the fields E and B associated with the accelerated 
particles, as calculated from Maxwell's equations. 
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In the 3-D simulations by ISiverskv"& Zharkoval (l2009h . the y-direction is chosen to be 
invariant. The system is periodic in the z-direction, so that a particle leaving the system 
through a boundary far from the reconnecting current sheet (see bottom of Figure [34l l reap- 
pears at the opposite boundary. The current sheet half-thickness is taken to be to 100 cm and 
the period along the x-direction is 2000 cm. In the PIC simulations discussed here, the com- 
putation region is extended well beyond the null point in order to explore the effect on the 
acceleration of different magnitudes of the transverse magnetic field. Plasma is continuously 
injected from the x = ±1000 cm boundaries of the simulation region. 

In order to avoid numerical instabilities in the PIC method, constraints on the steps in 
time and space must be satisfied (for details, see ISiverskv & Zharkoval '2009Y To satisfy 
these conditions while keeping the code running time reasonable, an artificially low plasma 
number density n = 10'" — lO'^ m^^ is used, and the proton-to-electron mass ratio is artifi- 
cially reduced to irip/me = 100. The spatial simulation grid has from 10 to 100 cells in the z 
direction and 100 cells in the x direction with Az = Ad (the Debye length) and Ax = Ad/5. 
There are ~100 particles per cell on average and the time step is 6 x 10^'" s. 

4.8.2 PIC simulation results 

The results of a simulation carried out with an extremely low number density n = lO'* cm^^ 
are shown in Figure l4~T9] for both protons (left panel) and electrons (right panel). This low- 
density simulation was performed mainly to verify that the PIC code can reproduce the 
results obtained in test-particle simulations (see Section l4~4l Figure l4~6l l. 

Consistent with the test-particle results, the formation of separate beams of accelerated 
protons and electrons is clearly observed, and electrons and protons are ejected into the 
opposite semispaces with respect to the x = midplane. The protons are largely unmagne- 
tized inside the current sheet, so that the ejection velocity of the "bounced" and "transit" 
protons are almost the same. On the other hand, the electrons are much more magnetized, 
so that bounced electrons cannot reach the midplane and thus gain much less energy than 
transit electrons. In this particular case, however, the polarization electric field caused 
by the charge separation is rather weak (~20 V m^'). It must be noted, however, that the 
electron skin depth for this simulation is about 5000 cm, which exceeds by a large factor 
the size (100 cm x 1000 cm) of the simulation region in both directions. The i on inertial 
length , which is also instrumental in controlling the electron/proton decoupling jMandt et al.l 
Il994bl) . is about 1000 cm, comparable to the size of the simulation region. Thus artificial 
decoupling of electrons and ion motions does result in this simulation. 

Simulations for a higher density, « = 10* cm^-^, are shown in Figure l4~20l for protons 
(left panel) and electrons (right panel). For this density the electron skin depth is reduced 
to ~500 cm, which is within the size of the simulation region in either direction. Thus the 
decoupling of electron and ion motions is not an artefact of the size of the region compared 
to the skin depth. Although the assumed density is several orders of magnitude lower than in 
the solar corona (often taken to be 10'° cm^^), already several differences between the test- 
particle and PIC approaches appear. The dynamics of the relatively unmagnetized protons 
are not substantially different: their trajectories are close to those obtained from the test- 
particle simulations (Figure 14.6b . They are still ejected mainly into one semispace {x < 
0) with respect to the midplane, and their acceleration rate coincides with the theoretical 
magnitude given in Section [33] 

In addition, however, a small number of protons is ejected to the x > semispace, 
where electrons are normally ejected in test-particle simulations (Figure [4.201 left panel). 
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Fig. 4.19 PIC simulation snapsliots of low-density (n = lO'" m^^) plasma particles on the x-V- phase plane 
for protons (left plot) and electrons (right plot). The protons and the electrons enter the RCS from both sides 
and are ejected to opposite semiplanes either as a beam (transit particles) or as a flow ( bounced particles); see 
the te xt for details. The cun'ent-sheet parameters are the same as in Section l4^ (from lSiverskv & Zharkoval 
l2009h . 
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Fig. 4.20 Snapshots of higher-density plasma particles (lO'- m^^) on the x-V- phase plane (PIC simulations) 
for protons (left plot) and electrons (right plot) entering the RCS from both sides. Protons are still ejected to 
the same semiplane, while the electrons form a cloud which experiences multiple ejections and returns to the 
acceleration region. As a result, the majority of the electrons ends up being ejected to the same semiplane as 
the protons, with only a small fraction of electrons being ejected into the o pposite semiplane (see the t ext for 
details). The current sheet parameters are the same as in Sectionl4^(from lSiverskv & Zharkovall2009l) . 



Figure \4?2T\ shows the fraction of the protons ejected into the x > semispace as a func- 
tion of Bxo and B^; the depend ence on Byo quahtatively coincides with that obtained by 
IZharkova & Gordovskv'J ( |2004) . For small Byo, both PIC and test-particle simulations show 
that the protons are ejected symmetrically with respect to the midplane. However, in the 
PIC simulation, all the protons are ejected into the one semispace (x < 0) whe n Byo > 7 G 
(w B;o), while in the test-particle simulation jZharkova & Gordovskvvl 2004) . the accel- 
erated particles are fully separated if Byo > 1.5 x IO^^B^q. Thus, in PIC simulations, the 
particle trajectories have less asymmetry than those in the test-particle approach. 



4.8.3 Polarization electric field induced by accelerated particles 

The PIC simulations show that the induced magnetic field B is typically much smaller than 
the background B. On the other hand, inclusion of the induced electric field E is essential. 
Indeed, its absolute value is much larger (by an order of magnitude) than the field Ey induced 
by the reconnection itself (see Sections [4.4l and l4.7.1b . 

Figure l4~22l shows the value of the charge separation electric field perpendicular to 
the current sheet, averaged over the z coordinate, for different values of B^q and Syo- The 
field becomes stronger when either B^q decreases or Syo increases. Also, when Syo 0, the 
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Fig. 4.21 The fraction of protons that are ejected into the x > semispace, Dashed line: fraction vs. Byo 
(B\0 = 0.2 G). Solid line: fraction vs. Byo (Byo = 1 G). All other parameters of the current sheet are the same 
as in Figure |4^2Q1 



4000 



3000- 




-3000 



Fig. 4.22 Electric field induced by particles in the PIC simulations for different values of B,o and B,o- All 
other parameters of the cun'ent sheet are the same as in Figure ITZOl 



system becomes symmetric, and accordingly so does E^. Th e distribution of the associate d 
charge density p(jc) = (1/4ti) dE^/dx has been presented bv lSiverskv & Zharkoval ( l2009h . 



4.8.4 Particle trajectories 



As discussed above, proton trajectories do not change much in the presence of the polariza- 
tion field, the only difference being that the "bounced" protons have somewhat larger orbits 
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dSiverskv & Zharkoval |200^. leading to two small er peaks at about ±1 m in the charge 
density plot (Figure 7 in ISiverskv & Zharkoval[200^ . 

The trajectories of electrons are much more complicated. First, electrons which enter 
from the jc < semispace have dynamics similar to those of the "transit" electrons: they 
drift towards the midplane, become accelerated and are ejected to the jc > semispace. 
However, the polarization field E^ix), which extends well beyond the current sheet thickness 
and has a component parallel to the magnetic field, acts to decelerate the ejected electrons. 
For the chosen magnitudes of and By, the majority of electrons cannot escape to the x > 
semispace; instead, they are dragged back multiple times to the current sheet and after some 
time become indistinguishable from the electrons that entered from the x> Q semispace. 

The electrons that come from the x>Q semispace demonstrate rather different dynamics 
compared with the case when = 0. These electrons, which are "bounced" from the RCS 
in the absence of an E^, can now reach its midplane. In the vicinity of the midplane, electrons 
become unma gnetized and oscillate with a gyro-frequency determined by the magnitude of 
By (Figure 8 in lSiverskv & Zharkoval2 0()9). After a large number of oscillations, the electron 
is ejected. If the electron's irutial velocity is small, it can become quasi-trapped inside the 
RCS, similar to the situation in the vicinity of an X-type null point jPetkaki & MacKinnonI 
Such electrons are accelerated on the midplane, ejected from it, then decelerated 
outside the RCS and return back to the midplane. This cycle is repeated until the electron 
finally gains enough energy to escape the RCS. Since the magnitude of the polarization 
field Ex^{x) is smaller at ;c < than at jc > (see Figure l4~22l i. it is easier for the electron to 
escape to the < semispace. Thus, eventually, most of the electrons are ejected to the same 
semispace as protons, contrary to the results of the test-particle simulations, while creating 
an electron cloud around an RCS before their joint ejection with protons. It seems that the 
asymmetry condition is shifted to much high er magnitudes of the guiding fie ld than those 
found in the test-particle approach (>0.01B20, IZharkova & Gordovskvvl[2004) . 



4.8.5 Energy spectra 

As shown bv ISiverskv & Zharkovd ( l2009l) , the ejected electrons, in general, form a wide, 
single-peaked energy distribution with the width being of the order of the mean energy. This 
is different from the test-particle simulations, in which two narrow-energy electron beams 
are formed. Consistent with Section l44l the mean electron energy still strongly depends on 
Sj, with lower corresponding to a higher mean energy and a wider energy distribution. 



4.8.6 Turbulent electric field 



Another interesting effect observed in the PIC simulations for the selected physical condi- 
tions and proton-to-ion mass ratio is the excitation of plasma waves. 

As can be seen from Figure l4~23l the E- component of the induced electric field is struc- 
tured with a characteristic length scale of Awave « 2 m. This structure propagates in time 
in the positive z-direction at a speed Vwave « 1.3 x 10^ m s^'. This corresponds to an os- 
cillation period r„ave « 1.5 x 10^^ s, consistent with the plasma frequency for a density 
n = 10^ cm^-'. Note, also, that the oscillating component of the excited wave E~ is parallel 
to the direction of propagation which corresponds to the polarization of the Langmuir wave. 
The instability mechanism is evident from study of the electron velocity y,. As shown in 
Figure 14.241 the electrons have an unstable "bump-on-tail" distribution. The range of ve- 
locities Vj; for which the derivative df/dV~ is positive is from 1.3 x 10'' to 2 x 10'' cm s^', 
corresponding to the phase velocity of the associated Langmuir wave Vwave- 
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Fig. 4.23 Electric field E. induced by particles in PIC simulation (B,o = lO^^T", B,o = lO^'^T, B^o = 4 ; 
lO^^r, Eyo = 250 V m-', nip/m^ = 10, « = 10* cm-^). 
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Fig. 4.24 Distribution function of electrons calculated for the region with turbulence plotted in Figure 14^23 



The use of plasma instabilities to produce anomalously high values of resistivity in 
the reconnection region have been studied by many authors. A comprehensive review 
of these results is beyond the scope of this article; however, a brief commentary is ap- 
propriate. Both 2-D and 2.5-D PIC simulations have revealed various types of instabili- 
ties to be present, fro m the lower hybrid drift instability (LHDI ) leading to the Kely in- 
Helmholtz instabiUty (Landau & Lifshitz"l960': 'Chen et al.''l997': 'Lap enta et al.ll2003h to 



kink (Pritchett et al. 1996; Zhu & Winglee 1996; Daughton 1999) and sausage teiichnej 
1 19991) instabilities. However, full 3-D PIC simulations (wit h an artificially red uced speed of 
light c = ISVac and proton-to-electron mass ratio of 277. 8; IZeiler et af]|2002l) . while retain- 
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ing the lower hybrid drift instability, reveal no evidence of kink or sausage instabilities. This 
indicates that the appearance of various instability modes may be dependent on the details 
of the simulation used, and hence argues for an abundance of caution at this time in applying 
these results to actual flares. 



4.9 Particle acceleration in collapsing magnetic islands 
4. 9. 1 Tearing-mode instability in current sheets 



As mentioned in Section [321 in reconnecting current sheets, diffusion converts magnetic 
energy into Ohmic heat, with a concomitant slow expansion of the current sheet at a rate 
proportional to the resistivity. Consequently, one way of enhancing the reconnection rate is 
to conside r enhancement of the resistivity through resistive instabilities such as the tearing 
instability jpurth et al.lll963h . 

In a current sheet wide enough to meet the condition TiHf 2> Ta, where Ta = cI/Va is the 
Alfven traverse time, the tearing instability occurs on a timescale of ldif{'^A/'^dif)^ ^ with the 
parameter A in the range < A < 1 . The tearing instability occ urs on wavelengths gr eater 
than the width of the current sheet, {ka < 1), with a growth rate jPriest & Forbe sl l2000h : 

® = (4.20) 

where the wave number k satisfies the condition [idif/'^A)^^'^ < ka < 1. The fastest growth 
is achieved for the longest wavelength (i.e., at ka = 1), while the smallest wavelength of the 
instability grows more slowly, on a timescale fj/j T^^^. As first demonstrated bv'F urth et al.l 
(|l963) (see also Spicer 1977), tearing-mode instabilities generated during reconnection in a 
current sheet can produce closed magnetic "islands" with an elliptical O-type topology. Ac- 
cordingly, large current sheets often observed in solar flares (e.g., see^ Sui & Holmanl2003l 



and references therein) can, depending on their sizes and conditions, be subject to the forma- 
tion of magnetic islan ds . For example, for a current sheet of size ~ 1 0,000 km (as reported by 
ISui & Holma"3l2003[) . embedded in a (coronal) medium with magnetic field S w 100 G and 
density n = lO'" cm^-', then the Alfven velocity is about 2000 km s^' and Ta w 5 sec. With a 
rather conservative collision frequency estimate 10^ s^', Tj,y w 10'^ sec. Substituting these 
values into Equation 14. 201 gives an island formation time in the range 10 — 1000 s. 

4.9.2 Studying particle acceleration in magnetic islands using a PIC approach 

In the past decade, th e tearing-mode in s tability has been extensively investigated, both in 

2.5-D current sheets dOrake et al.l ll997: Shav & Drake 1998; Drake et al. 2005, 2006b','a) 

and in a full 3-D kinetic PIC approach tZeiler et al. 2002; Pritchett & CoronitL2004. ; ,Karimabadi et al.l 
l2007l;IShav et alJl2007l;lD"aughton et al.ll2009l;lDrake et al .ll2010l) . 



lShav& Drake' ('1998') studied the dynamics of whistler waves generated inside the mag- 
netic islands, using a code in which the electric and magnetic fields were calculated at each 
time step using the Poisson condition V -E = 4Kp . The time step was normalized to the elec- 
tron cyclotron time Q^^, calculated for a maximum initial magnetic field Bq , the spatial step 
normalized to the skin depth c/cOpe, and the velocity normalized to the Alfven velocity for 
electrons Vac = Bq j \f^%n^me. The speed of light was artificially reduced to SVap in order to 
marginally resolve Debye radii. The simulation grid was 512x512 with 7 million particles 
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Fig. 4.25 Out-of-plane current density in 2.5-D PIC simulations of reconnection with a guiding field, for 
four consecutive time intervals, from top to bottom. The environment contains two current sheets and is 
characterized by the appearance of islands formed by the tearing instability. At the start, there are only (a) 
large islands, followed at later times (b-d) by the formation of smaller magnetic islands. Note that such islands 
do not appear in 2-D PIC simulations of the same environment without a guiding field. From iDrake et"^ 
j2006bh . 

involved. Small magnetic perturbations and associated currents were added as initial con- 
ditions to form seed magnetic island s centered around the two current layers, as shown in 
Figure |4~25] Using such a formalism. IShav & D rake ( 1998) were able to resolve the electron 
dissipation region, despite the use of an artificially low ion-to-electron mass ratio of 200. 

A further modification of the PIC code was used to extend it to a 3-D model with 64 
grid cells in the third direction, whil e increasing the p roton-to-electron mass ratio to 277.8 
and the speed of light to c = XiV^e dZeiler et al.ll2002h . Unlike a conventional Harris equi- 
librium, in which the variation of the density produces the equilibrium current, the density 
was taken to be nearly constant, with the current being produced by the E x B drift of 
the electrons. This resulted (see Figure l4~26l l in a "parallel" electric field occurring in the 
plane of the island current s heet iDrake et al. 1997: Shay & Dra ke 1998; Zeiler et al. 2002). 
IPritchett & Coronitil ( l2004h have further shown that in 3-D current sheets with a guiding 
magnetic field, the "parallel" electric field is formed only in the two quadrants {-\-x^ -\-y) and 
{—x^—y). Electrons entering these sectors are accelerated by this "parallel" electric field; 
the other quadrants are characterized by an excess of ions, strongly reminiscent of the drift 
separation of electrons and protons in a 3-D single current sheet (Sections I4.4l and l4~8] >. 

Once formed, magnetic islands tend to shrink due to magnetic tension forces. Conserva- 
tion of the longitudinal invariant leads to efficient first-order F ermi acceleration of particles 
tied to the island's field lines. As shown bv lDrakeetal] ( l2006cJ) , the process of island shrink- 
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Fig. 4.26 Current density in PIC simulations for times (a) t = 1 .1 / cOci 
correspond to magnetic flux lines. Adopted from lZeiler et al.l )2002l) . 



and (b) t = 12.9/CB„. The solid lines 




Fig. 4.27 The multi-scale structure of the electron diffusion region around an X-type null point of total length 
of Ag and thickness 5^, as derived from 3-D PIC simulations. The solid lines show the electron trajectories 
for the time ta>ci = 80 {cOd is the ion cyclotron frequency). In the inner region of length Dg, there is a steady 
inflow of electrons and a strong o ut-of-plane current; the outer region is characterized by electron outflow 
jets. From lKarimabadi et aljj2007l) . 



ing halts when the kinetic energy of the accelerated particles is sufficient to halt the further 
collapse due to magnetic tension forces; the electron energy gain is naturally a large frac- 
tion of the released magnetic energy. A sufficient number of electrons can be produced to 
account for a moderately-sized solar flare. The resultant energy spectra of electrons take the 
form of power laws, with harder spectra corresponding to lower values of the plasma j3. For 
the low-jS conditions appropriate to the solar corona, the model predicts accelerated spectra 
that are much harder (5 w 1.5) than observed. 

Further progress with 3-D PIC simulations of magnetic islands, using a larger 2560 x 
2560 grid, revealed that the electron diffusion region formed by the electron flows occurs 
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Fig. 4.28 Evolution of the cun'ent density over a region of extent 800^/, (where d, is the ion inertial length), 
for different times t. White hnes are the magnetic flux sur faces. The bottom plot presents a close-up of a 
region in which formation of new islands is occurring. From lDaughton et alj i2009h . 



over a width of 5de (where de is the electron inertial length -'Bim et al."200l'), while other 
simulations show an ex tent up to lOdj (where dj is the ion inertial length - Karimabadi et alJ 
l200llShav et"al]|2007l) : see Figure 14.271 The electron diffusion region becomes elongated 
with a two-scale structure: an inner region (size in Figure l4~27l i with a strong electron- 
drift-associated electric field (see Equation 13.911 and a wi d er diffusion region characterized 
by electron jets ( iKarimabadi et al.ll2007l : [shav et al.ll2007l) . IShav et al.l (120071) found that the 
electrons form super-Alfvenic outflow jets, which become decoupled from the magnetic 
field and extend to a large distance from the X-type null point, similar to the manner in 
which electron jets decouple from the protons in a test-particle approach to reconnection in 
a single 3-D current sheet (Section l4~4l i. 

Recently, lOauehton et al ] il2009l) showed that for coUisional Sweet-Parker reconnection 
with Lundquist numbers S higher than 1000, the reconnection region becomes susceptible to 
the tearing instability, leading to the formation of many magnetic islands with half-thickness 
approaching the ion inertial length di (see Figure l4~28l l. This result shows that the current 
sheet size is a key factor in defining timescales and the role of different physical processes 
at different stages of the reconnection process. For example, the shrinking of a large cur- 
rent sheet during MHD reconnection is halted when it approaches a size conducive to the 
formation of magnetic islands. 
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4.10 Limitations of the PIC approach 



Many PIC simulations (IjBim et al. l200lMDrake et alj2006bl : lTsiklauri & Harukil 20071 : 1 Shay et alj 



bOoilKarimabadi et alJl20o¥ have been concerned with the study of reconnection rates 



in 

the vicinity of a magnetic X-type null. Even for the very low densities used in these sim- 
ulations, it is apparent that the feedback of the fields induced by accelerated particles can 
be substantial (see Section |4~7] |: it can significantly modify the particle trajectories and the 
energies gained during their acceleration. 

However, because of computer power limitations, the size of the region in which the 
induced electric and magnetic field are simulated and coupled to the particle motion is neces- 
sarily rather small: ~ 100 x 50 x 1000 cm for 2-D simulations ( Bim et alj200ll:lTsiklauri & Harukil 



2007h. or 4000 x 40000 x 86 cm for the most recent 3-D simulations dShav et alJl2007l : 



Drake et alJ2"oioh . These are orders of magnitu de smaller th an the 10^ - 10^" cm sizes asso- 



ciated with flares ( iMasuda et al.lll994[Sui & Holman,2003l : Ijiang et al.ll2006) . Only in the 
PIC simulations discussed in Section |4~8] is the computation region extended well beyond 
the null point in order to explore the variations of acceleration process, at different magni- 
tudes of the transverse magnetic field. Thus, the narrowness of the simulation region still 
imposes strong limitations on the number of particles being accelerated in PIC simulations. 

Other puzzling features of the simulations relate to the origin of the electric fields in 
the models. One example is the "parallel" electric field reported in many 3-D simulations 
dZeiler et al.|[2002l : IPritchett & Coronitill2004l : iDrake et allHoioh and produced either as a 
result of an advective electric field (see the first term in Equation |3.9l l and/or by a gradient 
of the magnetic field component across the current sheet (see the second term in Equa- 
tion [3]9l This 'parallel' electric field is the electric field induced by accelerated electrons 
known as a self-induced (return-current) electric field associated with their propagation into 
a flaring atmosphere, which can exceed the local Dreicer field by a factor of up to ~30 
Izharkova & Gordovskvvll2006h . Although both the "parallel" electric field, or the "self- 
induced (return current')' electric field associated with the propagation of accelerated elec- 
trons can cause various interesting effects (e.g., the charge separation electric field can decel- 
erate electrons and even return them back to the acceleration region), causality arguments 
would strongly suggest that neither field can be considered as the main driver of particle 
acceleration since the field i& formed by the accelerated particles themselves. 

A related problem is the difference in electron and ion gyration radii/periods ( iBirn et al] 

l200lh . Because of this difference, within a timescale of < 10 ^ s, electrons can move to- 
wards the midplane, be accelerated and ejected, while the protons have only just started 
their gyrations in the RCS environment. In the simulations to date, the proton-to-electron 
mass ratio has been arbitrarily lowered to 100-200, an order of magnitude smaller than the 
actual value. Even the speed of light is conveniently reduced by some one or two orders o f 
magnitude to allow a reduction in simulation times ( iBim et al ] |200ll : brake et al.l l2006cJ) . 
Such simulations, then, serve to highlight that by artificially adjusting fundamental physical 
constants (such as the electron mass or the speed of light) in order to achieve results with 
a reasonable amount o f compu tational effort, certain features of the results can be altered 
significantly. Although lShav & Drake ( 1998) have shown that the reconnection rate may be 
indeed insensitive to the adopted values of both the electron mass and speed of light. 

Issues such as the ratio of macro- to micro-scales can be resolved only by increasing the 
reconnection region size in the PIC simulations up to sizes comparable with observations. 
However, even with exponential increases in computing power, there are substantial doubts 
that one can achieve this goal within the next few decades. Possibly a hybrid approach, com- 
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billing PIC simulations convoluted with test-particle extensions to greater source volumes, 
will allow us to make substantial progress in this area. 

5 How well do Recent Theories Account for Observed Flare Characteristics? 

5.1 Interrelation between acceleration and transport 

Before we address the relative merits of the various particle acceleration models, it is im- 
portant to realize that many flare observations consist of radiation integrated over the source 
volume; the quantity observed involves integration over not only space but also quantities 
such as particle pitch angles. Therefore, in order to relate to observed features (Section]!}, 
one must incorporate the effects of particle acceleration with the transport of the accelerated 
particles throughout the remainder of the flaring source. Although such an inclusion will in 
all probability obscure information on the acceleration process, it can also, more profitably, 
be used as a diagnostic of particle transport processes themselves. 

As we have seen (Sections 14.41 14.51 and 14.6b . accelerated electrons can attain quasi- 
relativistic energies within an extremely small timescale and can thus propagate over the 
extent of a solar flare of size ~10^ cm within '-^O.Ol-O.l s. The proton acceleration time is 
longer than that for electrons by up to two orders of magnitude (Section 14.4b . so that pro- 
tons still remain within the acceleration region while the electrons are ejected from it. Thus, 
particle transport mechanisms ar e essential for cornpletin g the electric circuit for particles 
escaping the acceleration region ( Emslie & Henoux 1995) and/or in establishing co-sp atial 
return currents ( Kni g ht & Sturrock 1977; Emslie & Henoux 1995; Zharkova & Gordovs kvv 

2006f). Such return curre nts can be formed both by ambient electrons ( Knisht & Stur rock 

i977':'Noc era etal.ll985h and by relatively high-energy beam electrons ( Zharkova & Gordovskvvl 
2006) that return to the corona with nearly the same velocities as the injected electrons. 
The self-induced electric fields created by the beam electrons and/or the return current 
have a structure very favorable for the generation of plasma turbulence ( iKarlickv & Bartal 
l2008h . which naturally creates the conditions required for stochastic a cceleration of elec- 
trons to sub -relativistic energies as discussed in Section 13.41 (see also lMilleretalj[l99^ ; 
IPrvadko & Petrosian.1997) . 

Another issue that can arise is consideration of radiative transport effects, such as ra- 
diative transfer of optically-thick lines and continua, and reflection (or albedo effects) of 
emission of hard X-rays from the photosphere. 

5.2 Testing acceleration models against observational constraints 

Here we briefly summarize the extent to which the various models discussed in Section |4] 
have increased our ability to account for the characteristics (e.g., number, spectrum) of ac- 
celerated particles in a solar flare. In other words, can any hitherto-identified acceleration 
process produce the rate, number and energy spectrum of accelerated particles inferred from 
observations? 

5.2.1 Stochastic acceleration 

Stochastic acceleration remains one of the most promising candidates to account for most of 
the characteristics deduced from hard X-ray, radio and gamma-ray observations. This mech- 
anism does, however, require a pre-supposed level of either plasma or MHD turbulence, the 
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basis for which has not yet been fully established. Further, high-frequency waves near the 
plasma frequency can be excluded as drivers for stochastic acceleration, since they would 
couple into a high level of decimetric radio emission in every flare, contrary to observations 
teenz et alfcOOSt) . 

On the other hand, acceleration by transit-time damping is compatible with a frequent 
lack of radio emission because the frequency of the postulated MHD turbulence is far be- 
low the plasma frequency. The non-thermal electron distribution "grows" out of the thermal 
population, resulting in a "soft-hard-soft" evolution of the hard X-ray spectrum that is con- 
sistent w ith observations. Further , such "soft-hard-soft" behavior has been found in coronal 
sources, teattaglia & Benj|2007h and is therefore inherent to the acceleration process and 
not to collisional transport effects. RHESSI observations showing the existence of a "pivot 
point" energy at which the differential electron flux does not vary with ti me, leading to the 
observed "soft-hard-soft" behavior, have been used ( iGrigis & Benzll2004h as evidence for a 
stochastic acceleration model. 

Recent models of stochastic acceleration that in volve the interaction of tu rbulence with 
high energy particle beams injected from elsewhere jSykov & Fleishmarill2009i) can also ac- 
count for some basic features of accelerated particles and their evolution with time. During 
the initial (linear) phase of the acceleration, effective particle acceleration by the longitudi- 
nal large scale turbulent motions, leading to visible spectral hardening. Then, because the 
accelerated particles eventually accumulate a considerable fraction of the turbulent energy, 
they deplete the energy in the turbulence, leading to a spectrum softening. Again, this is 
consistent with observations (Battaglia & Benz 2007) of "soft-hard-soft" spectral evolution 
in the main stage and even some hardening (soft-hard-harder behavior) during later stages 
jKmcker et al.ll2007ah . 

5.2.2 Collapsing current sheets 

Estimates of the electron energy gained in collapsing magnetic traps associated with recon- 
nection in a helmet-type geometry, with and without a termination shock, reveal that there 
is a sufficient amount of energy released in accelerated particles during such a process. In 
addition, betatron acceleration, possibly coupled with first-order Fermi acceleration in col- 
lapsing current sheets, can account for the electron number and energies gained, although 
the energy spectrum of the accelerated particles resembles more the higher energy tail of a 
Maxwellian distribution. It remains to be established how the accelerated particles can move 
across the magnetic field of the trap in order to precipitate into the flare footpoints. 

5.2.3 Acceleration in 3D reconnecting current sheets 

New MHD simulations of particle acceleration rates and spectra from collapsing current 
sheets have revised our views on the efficiency of particle acceleration inside a diffusion 
region, although the number of particles (both electrons and protons) that can be accelerated 
in a current sheet i s limited by the volume encompassed by the current sheet ( iHolmanll 19851 : 
iLitvinenko' 1996a, see Section [331 l. However, more complex 3-D current sheets and/or mul- 
tiple current sheets may be reconcilable with this model (see Sections [4.4l and l43T >. Acceler- 
ation of electrons in a 3-D RCS occurs mostly around the RCS midplane, over a very short 
timescale comparable with the growth rate of hard X-ray emission in solar flares. Accelera- 
tion of protons has a much longer timescale and results in a spread of protons in a location 
outside the midplane. The particle orbits in the RCS lead to a separation of electrons and 
protons into the opposite semiplanes, corresponding to ejection in opposite directions and 
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into opposite loop legs. The separation of electrons and protons during acceleration in an 
RCS leads to the formation of a polarization electric field (induced by Hall currents), which 
sets the frame for the generation of strong plasma turbulence that can lead to further particle 
acceleration by stochastic processes. Single 3-D current sheets with guiding fields have been 
found, through both test-particle and PIC approaches, to be able to supply sufficient energy 
into both electrons and protons in order to account for RHESSI hard X-ray and gamma-ray 
observations. 

Spine reconnection in 3-D MHD reconnection simulations can naturally produce two 
symmetric jets of energetic particles escaping along the spine, while the fan reconnection 
produces azimuthally localized ribbons of high energy particles in the fan plane. Both fea- 
tures, in principle, are often observed in flares (e.g., magnetic jets in SOL2000-07- 14T10:24, 
X5.7, and and SOL2002-07-23T00:35, X4.8) and two-ribbon flares as discussed in lFletcher etaP 
( I2OIII) . 

Another possibility involves tearing instabilities in current sheets, which leads to the 
formation of large numbers of magnetic islands, on which a large number of particles can 
be rapidly accelerated; however, the simulated electron energy spectra are much harder than 
inferred from hard X-ray observations. 

Despite these qualified successes, there are two outstanding problems with this kind of 
acceleration model: (1) incorporation of the self-consistent electrodynamic fields associated 
with the accelerated particles; and (2) production of sufficient numbers of accelerated parti- 
cles to match the observations. The latter can be marginally accounted for if one notes that 
strongly stressed 3-D current sheets, discussed in Sections [4. 4l and l431 are characterized by 
much larger diffusion regions than in simple reconnection models (Section [T2l ). 

5.2.4 Complex reconnection models 

The other way to tackle the particle number problem involves multiple dissipation sites asso- 
ciated with self-organized criticality (SOC) that could lead to "spreading" of the dissipation 
across a large coronal volume (see Section l431 l. Particles move between the current sheets, 
either gaining or losing energy in each one, resulting in a systematic second-order gain in en- 
ergy. Particle energy distributions from the simulations fit the observed spectra for electrons 
very well, although proton spectra have yet to be simulated. 
However, a few essential issues do remain: 

1. The models invoke an ad hoc distribution of RCS sizes and physical conditions in order 
to account for the number and spectra of the accelerated particles; 

2. The basic model of a current sheet (discussed in Section [331 l is used to derive the energy 
gains by accelerated particles in a component current sheet, while the actual gains are 
significantly different in more realistic 3-D models of reconnecting current sheets; 

3. Even if multiple current sheets are formed in the whole volume of the corona, there 
remains the issue of where the reconnecting magnetic field lines are embedded and how 
the test particles can move across the strong magnetic field threading the reconnection 
volume in order to reach another current sheet; 

4. Some of these models also suffer from problems associated with the neglect of the feed- 
back of the electric and magnetic fields associated with the accelerated particles them- 
selves; 

5. Multiple RCS models simulated with a more realistic PIC approach are applicable to 
only a very small region (about 10^ cm^), and it is assumed without justification that the 
results can be scaled to flare volumes. 
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6 Summary 

Over the past decade, and especially in the RHESSI era, significant progress has been made 
in advancing various particle acceleration models and assessing their ability to account for 
the observed features in solar flares. RHESSI observations have been instrumental in encour- 
aging the further development and exploration of both "traditional" and "non-traditional" 
mechanisms, including more complicated magnetic field geometries and topologies, to en- 
hanced MHD and PIC models of magnetic reconnection in a single current sheet, to multi- 
ple current sheets and stochastic acceleration by various plasma waves. Substantial progress 
has been made through the use of phenomenological estimations made in collapsing traps, 
termination shocks of Alfven waves, and direct simulations using both test-particle and self- 
consistent PIC approaches. 

Consideration of the full three-dimensionality of the magnetic reconnection process has 
been found to significantly enhance the acceleration efficiency in current-sheet diffusion re- 
gions. However, the question still remains whether such a mechanism can account for the 
very large number of electrons responsible for hard X-ray emission and whether consider- 
ation of the whole electric circuit formed by precipitating and returning beams in recon- 
necting loops can resolve the particle number problem. PIC simulations have demonstrated 
that the acceleration of substantial numbers of charged particles causes a number of im- 
portant secondary effects: for example, polarization electric fields induced by separation of 
electrons and protons towards the midplane, return currents formed from precipitating and 
ambient electrons, and plasma turbulence related to beam instabilities. 

Interactions of large-scale loops with strong magnetic fields can lead either to accelera- 
tion in collapsing traps by betatron and/or first-order Fermi mechanisms, or to acceleration 
in the termination shock formed by such the collapsing sheets (or by large-scale Alfven 
waves produced during the reconnection). Interaction of smaller-scale magnetic fields can 
produce magnetic islands or multiple current sheets in which particles can be also accel- 
erated in a cascade processes if the particle travels and gains energy from one feature or 
another. 

The exact acceleration scenario will depend on the magnetic field topology and on the 
absolute magnitudes of physical parameters. Further, it is possible, or even likely, that all 
the acceleration processes discussed above play some role in solar flares and elsewhere in 
astrophysics. The strict demands imposed by the RHESSI observations have necessitated a 
thorough and critical review of all these models, leading to substantial new understanding 
along the way. 
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